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Abstract: We analyze the large frequency behavior of the spectral densities that 
govern the generalized Langevin diffusion process for a heavy quark in the con- 
text of the gauge/gravity duality. The bare Langevin correlators obtained from the 
trailing string solution have a singular short-distance behavior. We argue that the 
proper dressed spectral functions are obtained by subtracting the zero-temperature 
correlators. The dressed spectral functions have a sufficiently fast fall-off at large 
frequency so that the Langevin process is well defined and the dispersion relations 
are satisfied. We identify the cases in which the subtraction does not modify the as- 
sociated low-frequency transport coefficients. These include conformal theories and 
the non- conformal, non-confining models. We provide several analytic and numerical 
examples in conformal and non-conformal holographic backgrounds. 
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1. Introduction and results 

Heavy-ion collision experiments and related data on the deconfined phase of QCD, 
[0] , have provided a first window for string theory techniques to meet the real world. 
The context is strong coupling dynamics near and above the deconfining transition in 
QCD. String theory via the AdS /CFT correspondence has provided a framework in 
order to understand strong coupling dynamics in the deconfined phase including the 
calculation of transport coefficients. Recent reviews on the progress in this direction 
are 0, |, |, §, §. 

Observables of particular importance are associated to heavy quarks. Heavy 
quarks may be produced in the Quark-Gluon Plasma (QGP) of the RHIC fireball 
and are then traveling to the detectors while moving through the dense QGP. They 
can be tagged reasonably well and are therefore valuable probes of the dynamics in 
the plasma and in particular for the mechanism of energy loss. 

A single heavy quark can be modeled in string theory by an open string. Its 
end-point is representing the heavy quark while the string is trailing behind as the 
quark moves. The large mass limit is important in order to neglect the non-trivial 
fiavor dynamics associated with light quarks (although with improved techniques the 
light quarks may also eventually be addressed reliably in the holographic context). 
As quarks are associated with strings ending on fiavor branes, a heavy quark ends on 
a brane that is stretching in the UV part of the bulk geometry. The motion of such a 
string, and the associated force acting on the quark from the thermal medium, have 
been studied in detail with several complementary methods, 0, |], |[ . In the simplest 
setup, the UV endpoint of a fundamental string is forced to move with constant 
velocity v along a spacial direction. The equations of motion for the full string are 
solved and the radial profile of the trailing string is found as it moves in a bulk black 
hole background representing the deconfined heat bath. The energy absorbed by the 
string is calculated and the drag force of the string is obtained. The picture remains 
roughly valid, while details change when conformal invariance is broken, 

An important improvement in this picture consists of the study of the stochastic 
nature of this system, in analogy with the dynamics of heavy particles in a heat 
bath giving rise to Brownian motion. This involves a diffusive process, that was ffrst 
considered in a holographic setting in by using the Schwinger-Keldysh formalism 



adapted to AdS /CFT in ^ 



Subsequently, a study of the (quantum) ffuctuations of the trailing string, [|I^, 
T3|, provided the information on the momentum broadening of a heavy quark as it 
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moves in the plasma. The stochastic motion was formulated as a Langevin process, 
0, , associated with the correlators of the fluctuations of the string. 

Many heavy quarks in experiments are relativistic. Therefore it is necessary 
to study the associated relativistic Langevin evolution of the traihng string, a feat 
accomplished in the Af = 4 case in and in general non-conformal plasmas in 
[ pOj , pT| . The same type of Langevin process was studied in for the case of an 
accelerating quark in the vacuum (rather than in a deconfined plasma), by analyzing 
the fluctuations of a trailing string in AdS with a non-uniformly moving endpoint. 

On the experimental front, there have been several results from the RHIC exper- 
iments, p2|-p6l. The experimental signatures are currently summarized by the 
spectra that originate in the semileptonic decays of charmed and bottom hadrons. 
From these spectra a modification factor i?^^ and an elliptic flow coefficient 
extracted. They capture the effects of the medium to the propagation of the heavy 
quarks. The data exhibit a substantial elliptic flow, up to ~ 10%, and a high-p^ 
suppression down to -R^^ ~ 0.25. These values are comparable to light hadrons. 
Radiative energy loss models based on pQCD, |^, do not seem to explain well the 
experimental data, . Elastic scattering energy loss plus non-perturbative interac- 
tions can on the other hand accommodate the data, |^ . 



In particular, the Langevin approach has been applied to the study of the heavy 
quark energy loss by several groups, and the related physics is summarized in the re- 
view . The Langevin evolution used was relativistic and with symmetric diffusion 
coefficients. As there was no microscopic model to provide the proper fluctuation- 
dissipation relation, the Einstein equations used vary, and in all examples it was 
assumed that the equilibrium distribution is the Jiittner-Boltzmann distribution. 
Moreover various combinations of friction forces were used, resonance models, pQCD, 
A^=4 AdS /CFT and combinations. A further recent analysis was performed in 



with similar conclusions. The associated relativistic and isotropic Langevin systems 
used have been introduced in the mathematical physics literature rather recently, 
[0 (see 1^ for a review). 

Before proceeding to the issues studied in this paper, we present briefly the 
results of |^ in order to put our conclusions here in context. In ||21|] a large class 
of non-conformal backgrounds captured by Einstein-dilaton gravity with a dilaton 
potential in 5 dimensions were considered in the context of heavy-quark energy loss. 
In a series of recent works, such backgrounds were analyzed both qualitatively and 
quantitatively and have provided a rich variety of holographic bulk dynamics. In 
particular, for a selected class of scalar potentials, they mimic the behavior of large- 
Yang Mills, [^-|^. This match can be quantitative, [^, agreeing very well both 



at zero and finite temperature with recent high-precision lattice data, |40]. However 



the Langevin analysis was done for generic non-conformal backgrounds in ^1 



The main method is to consider a string end-point that is forced to move with 
velocity v. Solving the Nambu-Goto equations of motion, the classical profile of the 



3 



trailing string can be found. The string stretches inside the bulk until it becomes 
completely horizontal at some value of radial coordinate r^, given by f{rs) = v"^ 
where /(r) is the blackness function of the background. When the quark is moving 
slowly, as f — )■ 0, the point approaches the bulk black hole horizon. 

The induced metric on the string world-sheet has the form of a two-dimensional 
black hole metric with a horizon at r = as first observed in |15[.^ This black 



hole is an important ingredient of the dynamics of the system. In particular it is 
crucial in the calculation of the thermal correlators using the Schwinger-Keldysh 
formalism, as well as for the fluctuation-dissipation relation. The world-sheet black 
hole has an associated Hawking temperature Ts that depends on several parameters: 
the background temperature T, the zero-temperature bulk scale A and the quark 
velocity v. It coincides with the temperature T of the heat bath only in the non- 
relativistic limit. In the conformal case, one has = Tgconf = ^(1 ~ v"^)^ < T. 
The numerical analysis of many examples shows that Tg < Ts^conf < T. The equality 
Ts = Tg^conf, in the first relation is attained, for arbitrary v, in the high T limit and 
also in the ultra- relativistic limit, f — )■ 1. 

The small fluctuations around the classical string profile satisfy second-order 
radial equations that are related to the associated thermal correlators by the holo- 
graphic prescription. It should be emphasized that such correlators are thermal with 
temperature Tg and not the temperature T of the heat bath. Moreover, they sat- 
isfy the fluctuation-dissipation relation associated with the emergent temperature 
Tg. The fact that the string fluctuations see a modified temperature crucially affects 
the Einstein relation between the diffusion constants. 

At the quadratic level of fluctuations, a relativistic Langevin diffusion equa- 
tion is obtained using the AdS /CFT prescription. The properties of this relativistic 
Langevin evolution differ substantially from rotationally invariant equations that 
have been introduced recently in mathematical physics [^, TL, In particular. 



here the evolution is not rotationally symmetric, and the Einstein relation is different, 
because the fluctuation-dissipation relation is different. This implies that the equi- 
librium configuration is not the standard rotationally-invariant Jiittner-Boltzmann 
distribution. 

The main focus of was on the local Langevin equation, which arises when 



looking at the large-time limit of the fluctuations of the heat bath, i.e. at the small 
frequency modes. On the other hand, the holographic computation gives access 
to the full frequency spectrum of the correlation functions driving the generalized 
Langevin dynamics. The holographic computation of the full Langevin correlators 
and the associated spectral densities for improved Holographic QCD was done mostly 



^This is a generic effect on strings and D-branes embedded in black hole/black-brane back- 
grounds. It was first observed in where it was used to propose that a different speed of hght 
is relevant for such branes. It is implicit or explicit in many holographic computations using probe 
flavor branes, p2, and strings p5|. 
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numerically in Analytic expressions (in terms of the bulk metric and dilaton 
profiles) were obtained in the two opposite regimes of small and large frequencies 
u (compared to an appropriate temperature scale) for general non-conformal back- 
grounds. 

In the large-frequency regime, the spectral densities are obtained via a modified 
WKB method, similar to the one followed in |^ for bulk fluctuations in an AdS- 
Schwarzschild background. The high-frequency behavior is different, depending on 
the mass of the probe quark. For finite mass, and for large u, the spectral densities 
grow linearly with u, whereas in the limit when the quark mass becomes infinite this 
behavior changes to a cubic power-law. However as we will see later in this work, a 
part of this regime is unreliable. 

Going beyond the zero-frequency limit is necessary when the diffusion process 
happens on time scales comparable to, or smaller than the auto-correlation time of 
the fluctuation propagators. More specifically, since these are thermal correlators 
at the temperature Tg, the large-time approximation breaks down over time-scales 
shorter than T~^. This condition puts a temperature-dependent upper bound on the 
momenta of the heavy quark, above which the diffusion process cannot be described 
by a simple local Langevin equation with white noise. In this context, it is useful 
to have approximate expressions valid for large frequencies (for examples, those we 
obtain with the WKB method) to model the behavior of the system in the regime 
where the local Langevin approximation breaks down and the dynamics becomes 
non-markovian (due to a non-trivial memory kernel) 

The numerical evaluation of the diffusion constants may lead directly to a com- 
parison of the jet-quenching parameters between the holographic QCD model and 
data. It was found that q± displayed a mild momentum dependence for large quark 
momenta, which however differs from the one obtained holographically in the con- 
formal case. As the temperature rises, q± increases significantly, approximately as 
~ T^. Interestingly, it was found that for temperatures above ~ 400 MeV, the lo- 
cal description of the diffusive process breaks down for charm quarks with momenta 
above ~ 5 — 10 GeV. This is because the process occurs on time scales shorter than 
l/Tg. This implies that in order to describe heavy charm quark diffusion in the AL- 
ICE experiment, one would need the full generalized non-local Langevin equation, 
and the full frequency-dependent correlator, rather than just its low-frequency limit 
captured by q±. This would constitute an interesting testing ground for holographic 
models, where the full correlators can be easily computed. The energies at which the 
energy loss mechanism described, is not any more the dominant one and radiation 
becomes the dominant mechanism were also determined. This was estimated by re- 



quiring that r<j remains below the would-be position of the flavor brane, ||15| , |5^ . It 
was shown that these limit do not substantially constrain this framework. 

The issues analyzed in this paper arose from the intent to apply the results of 
to comparison with experimental data. It was argued in that for parts of 
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the phase space of heavy quarks in LHC, the full non-Markov Langevin is relevant, 
and therefore the full Langevin correlators are needed. When trying to apply this 
Langevin evolution to real quarks, two immediate questions arise: 

• What is the proper behavior of the correlators at high-frequency u, and how 
is this expected to affect heavy quark diffusion? 

• What is the proper observable motion of heavy quarks in a strongly coupled 
plasma? 

A direct derivation of the Langevin evolution in the vacuum of the strongly 
coupled theory (at T = 0), involves the zero temperature retarded and symmetric 
correlators. As usual the retarded correlator controls the traditional dissipative term 
while the symmetric correlator the fluctuations of the Langevin noise. 

As shown in earlier works, |T^, |T^, the retarded T = correlator in AdS, 
is Gr ~ (in the infinite quark mass limit). The same is the large frequency 
behavior in any asymptotically AdS case, ||2l]]. This short distance behavior provides 
a Langevin evolution at T = with a "dissipative force" that contains the third time 
derivative of the coordinate. Obviously this behavior is unphysical and the short 
distance behavior of the Langevin evolution must be redefined so that they is no 
semiclassical dissipation in the ground state {AdS vacuum). 

A related issue concerns the causality of the evolution ingrained into the causality- 
related dispersion relations (|2.8| ) for the force correlators. Such dispersion relations 
are valid only if the spectral densities vanish sufficiently fast at large frequencies. 

The above considerations suggest that the proper definition of the Langevin 
diffusion of the heavy quark in the plasma must be defined so that the diffusion 
in the vacuum is trivial. This involves the redefinition (dressing) of heavy quark 
coordinates together with an appropriate modification of the path integral measure 
controlling their fluctuations. 

This procedure is developed from first principles in this paper. We find that in 
the regime in which the fluctuations of the quark coordinates are perturbative, the 
retarded correlator Gr{ijJ,T) is replaced with Guiuj.T) — Gr{uj,Q) which vanishes 
at T = and falls off (as 1/w) at large u. The corrections to this simple result are 
controlled by <^ 1. This condition is necessary because, in the gravity picture, 
having 70; ^ Mq implies a breakdown of the boundary effective description of the 
string endpoint as a free heavy quark. Also, in this regime non-linear backreaction 
effects become important. 

Similarly, the symmetric correlator is replaced by Gs{oJ,T) — Gs{u},0). the Ein- 
stein relation relating Gr to Gg gets modified, but it unchanged in the zero-frequency 
limit. Therefore, the relation between the friction and noise term in the local limit 
of the Langevin process is the same as in [PT[| . 
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A subtraction similar to the one we perform in this work has been recently 
advocated in the case of bulk correlators, in |^5|], in order for the corresponding 



spectral densities to satisfy appropriate sum rules. Here we are able to justify why 
this subtraction is physically sensible, in the case of the Langevin process, but we 
expect that similar arguments can be used for the bulk spectral functions. 

The need to define dressed correlators stems essentially from the need to cure a 
UV behavior which is too singular. On the other hand, we should expect that this 
does not affect the low-frequency regime, which gives the transport coefficients and 
should be independent of the details of how we treat short- distance physics. In other 
words, the subtracted correlators should have the same small frequency limit as the 
unsubtracted ones. 

In this paper therefore we perform the following task 



1. We define the dressed holographic Langevin evolution as described above, by 
performing the appropriate redefinitions in the path integral, and we show that 
they have the correct fall off at high frequency for the dispersion relations to 
hold. 

2. We analyze the low-frequency asymptotics of the zero temperature retarded 
force correlators in order to assess to what extent the dressing affects IR asymp- 
totics of the diffusion. 

3. We compute numerically the dressed correlators Gi{uj,T) — Gi{uj,Qi) both in 
AdS, and in exact scaling backgrounds as particular examples, and verify that 
the dressed correlator has the desired properties. We check explcitly that the 
dressed correlators have compact support and exhibit a few oscillations before 
they die-off beyond the natural scale set by Tg. 



In this paper we analyze the IR asymptotics of the retarded correlator of the force 
exerted on the fundamental string by the plasma as a function of the IR asymptotics 
of the dilaton potential. While the IR asymptotics are irrelevant for the large fre- 
quency behavior of the spectral densities, they are the feature that determines the 
diffusion constants. In this work we limit the analysis to models that give rise to 
zero-temperature geometries that are non- confining. In this case, the identification 
of the trailing string solution at zero-temperature is straightforward, and the dress- 
ing procedure can be defined consistently. The case of confining zero-temperature 
theories is of course more interesting , since it contains QCD, but it is also more 
subtle, and will be analyzed separately, [ |65[] . 

The structure of this paper is as follows. 

In Section ^ we review the Langevin dynamics of a relativistic point particle 
and its holographic description. We give general expressions for the large frequency 
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behavior of the spectral densities, and we discuss some exphcit examples. We extend 



the results found in |2T[ to non-asymptotically AdS scaling backgrounds. 

In Section ^ we compute the vacuum Langevin correlators, both in the high 
and low-frequency behavior, calculate the diffusion coefficients in various types of IR 
geometries, and give a classification of the IR geometries accordingly. 

In Section ^ we define the subtraction procedure from a fundamental perspective, 
using a path integral formulation of the Langevin dynamics. We show that the 
dressed correlators defined in this way arise naturally if one imposes some basic 
physical requirements. In this Section we also show that they have the correct fall- 
off behavior at large frequency to allow a consistent application of the dispersion 
relations. 

In Section 5 we present numerical results for the dressed correlators in both AdS 
and the non-conformal scaling geometries. 

Some technical details are left to the Appendices. 

2. Review of the Langevin evolution in strongly coupled plas- 
mas 

In this Section we review the description of generalized Langevin dynamics in holo- 



graphic finite-temperature plasmas, and the corresponding result |2l[ for the high- 
frequency behavior of the spectral density. These results are extended in the last 
Subsection to non-asymptotically AdS scaling backgrounds. 

2.1 Generalized Langevin dynamics of a relativistic particle 

A probe heavy quark propagating in a deconfined plasma undergoes a generalised 
Langevin process, due to the interaction with the medium. If we assume the quark 
to be heavy (compared to the temperature), we can instantaneously describe the 
fluctuations 6X^{t) of the quark position with respect to a straight trajectory, by an 
action of the form: 

S[SX{t)] = So[6X{t)] + j dr 5X^{T)J'^\r) (2.1) 

where 5*0 is the free particle action and J^^{t) is a field (operator in the operator 
formalism) which depends only on the medium degrees of freedom and dynamics, 
and represents the effect of the microscopic interactions with the plasma. Then, the 
generalised Langevin process is given by the equation: 

/ + 00 
rfrGt(r)5X,(t-r) + C{t). (2.2) 
-oa 

On the left hand side, P^{5X) represents the classical equation of motion for 5X^{t) 
in the absence of dissipation, derived from 5*0 In the absence of external forces. 
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P*(5X) typically takes the form of a two-time-derivative term: 

P\6X) = M'J^j6X„ 



(2.3) 



where ^^e// effective mass matrix. For a non-relativistic particle we simply 

have A^e// ~ MqS^^ , whereas in the case of a relativistic particle with velocity v, it 
is given by (see e.g. pl[]): 



Mr 



V V' 



^{v) = {l-v^Y^'^. (2.4) 



Mq being the quark mass. 

The effective force on the right hand side of ( |2.2|) is composed of two terms: 



an overall friction term, described by the convolution with a memory kernel 
given by the retarded correlator of the force operator J^{t): 



Gi{t)^-ie{t){[r{t),:FKm 



(2.5) 



a stochastic gaussian external force with moments given by the symmetrized 
correlator of the same operator J-": 



{m) = 0, 

(e(t)e(o)) = G%^{t) ^ {{r{t\7\^)}) . 



(2.6) 



It is important to stress that the expectation values are calculated in the appropriate 
ensemble. In the non-relativistic limit, this is the plasma ensemble (heat bath). 
However at relativistic speeds the ensemble is different from that of the plasma 
ensemble, and is characterised by an emerging temperature T, different from that of 



the plasma, [18, 21 



Going to Fourier space, one defines the spectral density by 



p^\u)= / dte'^' {[r it), T^iO)]). 

J — oo 

One then has the dispersion relations: 



(2.7) 



ReG%{co) 



P 



:2.8) 



Note also here that these dispersion relations are valid only for appropriately sub- 
tracted densities, and we will expand on this below. 

The retarded correlator does not depend of the specific density matrix of the 
medium, but only on the canonical commmutation relations (or OPE) of the opera- 
tors J-" (see e.g. EUI for a recent discussion). On the other hand, by specifying the 
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density matrix one obtains a relationship between the retarded and symmetric cor- 
relators. Thus, the spectral density plus a choice of ensemble determine completely 
the generalized Langevin equation. For example, in the case of a thermal ensemble 
in equilibrium at a temperature T, one has: 

Gsym = VT coth(a;/2T) p(w). (2.9) 

The crucial consideration that motivates the present paper is that, given a spec- 
tral density p{uj), in order for the real-time correlators to be well defined through 
equation (PT7|), and for the dispersion relations to be valid, p{u) must have a 
sufficiently fast fall-off at large u. If this is not the case, one has to perform a suitable 
subtraction in order to have a well defined Langevin process. As usual, a bad large 
frequency behavior is related to a too singular behavior of the real-time correlators 
in the coincident time limit. 

That this is precisely the case in our situation can be seen heuristically by di- 
mensional analysis: from equation ( |2.1| ), the force operator J^{t) has mass dimension 
2, therefore one expects that for short time and large frequencies one has, in a CFT: 

p{uj) ~ (2.10) 

In coordinate space this implies a correlator with a delta-function singularity at the 
origin, ImGRit) ~ 6"'{t). 

The holographic calculation shows indeed that the unsubtracted spectral density 
grows in fact like in the Mq — > oo limit, as we will review shortly^. With this 
behavior, the Fourier integral ( |2.7| ) and the dispersion relations ( p.8| ) are ill defined. 

The short-time divergence is determined by the UV behavior, which may be 
completely distinct from the dynamics of the medium in which the quark propagate. 
In the following sections we will define a subtracted spectral density, which physically 
means that the Langevin dynamics will describe a "dressed" quark, in which the effect 
of the vacuum contribution has been renormalized away. 

2.2 Langevin correlators in holographic plasmas 

In gauge/gravity duality, a heavy external quark moving through the plasma at tem- 
perature T can be described by a string whose endpoint at the boundary follows the 
quark's trajectory 0-||3^. The string extends into the bulk, whose geometry is a 



black hole background with appropriate temperature T. The bulk fields 6X^{r,t), 
describing the fiuctuation of the trailing string around its classical solution, corre- 
spond holographically to the force operator on the boundary. Thus, we can extract 



^For finite Mq, it grows instead linearly in uj, for uj ^ Mq. However, as we will argue below, 
this regime is unphysical, since in this case the holographic picture of a heavy quark as the endpoint 
of the string with boundary is inconsistent. Therefore, by "large frequency" we will always mean 
a; ^ T, and for finite quark mass, one should always restrict the analysis to T ^ a; ^ Mq. 
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the force correlators of the previous section from the behavior of the traihng string 
fluctuations, using the standard holographic rules. 

Here, we will consider five- dimensional planar black holes, with metric'^: 

f{r)dt^ + dx'dxi 



ds' = h\r) 



(2.11) 



S 



The trailing string dual to the quark is governed by the Nambu-Goto action, 
1 



NG 



d cry — det g, 



al3 , 



gai3 = g^.^daX'^dfsX' 



/i, z/ = . . . 5 
a,P = 0,1 

(2.12) 

where g^j^i, are the components of the bulk metric ( |2.11|) . We choose the gauge ^° = 
t, = r for the world-sheet coordinates, and consider a classical string embedding 
of the form 

X(t,r) = ivt + ar))^, (2.13) 

i.e. such that the endpoint r = moves at constant velocity v, and the rest of 
the string trails along in the bulk. The induced 2d metric on the world-sheet has a 
horizon at r = r^, defined by 



firs) 

with associated Hawking temperature T^: 



1 

4:71 



firs) firs 



46'(r,) ^ /'(r,) 



(2.14) 



(2.15) 



b{rs) ' fivs 

The fields dual to the force operator J^it) are the fluctuations around the trail- 
ing string solution |2.13| . It is necessary to distinguish fluctuations longitudinal and 
transverse to the background velocity v. The fluctuating trailing string is described 
by the embedding 



X(t,r) = {vt + ^{r) + SX'>^{t,r)) - + 6X^{t,r), v ■ 5X^ = 



(2.16) 



where ^(r) is determined by the classical NG action. Expanding the Nambu-Goto 
action to second order we obtain the quadratic action governing the fluctuations. In 
a convenient coordinate system, where the world-sheet metric is diagonal this takes 
the simple form: 



^2 = -: 



drdr 



i=l 



(2.17) 



^This is in the string frame. In the whole of this paper the metric wih always be in the string 
frame. For a conformal theory this does not make a difference as the dilaton is constant. But for 
the non-conformal theories we will consider where the running scalar is the dilaton, it makes a big 
difference. 
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where the kinetic operators are defined by 



with 



'f - ssf^"^. - (2.18) 







C = vb\r,), Z = h'^^0^^. (2.20) 

For a harmonic ansatz of the form 5X*(r, r) = e*'^'^5X*(r, cj), the equations fol- 
lowing from the action ( p.lTj ) are: 

dr [R dr {5X^)] + ^ 5X^ = 0, (2.21) 



dr 

where 



i^i? 9. (5X11) 



+ ^5^" = 0, (2.22) 



i? = v/(/-^2)(&V-^^2)- (2.23) 



The holographic prescription for the retarded correlator, computed with the 
diagonal induced metric ( p.l9[ ), is given by: 

Gniu) = - [^Ur, u;)g^'d^^n{r, u;)l,^,,^^^ • (2.24) 

Here \E'ij(r, w) denotes collectively the fiuctuations ^X", ^X-*-, solutions of equations 
( P.2I - O2D , and the factor Q^^ is the appropriate one from equation ( 2.18 ). The so- 
lutions ^i^(r, w) must obey the appropriate boundary conditions: unit normalization 
at the UV boundary, and in-falling at the world-sheet horizon, i.e 

^R(r,w) ~ ^/,(w) (r, -r)"3^% r ~ r^. (2.25) 

The UV boundary is taken to be r = in case we consider the quark mass 
infinite, or r = tq > for a finite quark mass Mq. The value rg is determined by 
calculating the mass as the Nambu-Goto action of a straight string extending from 
tq down to the deep IR region. For a heavy quark, rg is close to the boundary and 
Mq ~ l/rg. 

The spectral density is related by equation ( p.8| ) to the imaginary part of equation 
( |2.24 ). Since the latter is a conserved fiux, one can evaluate it at the horizon r<. rather 
than at the boundary. The result is: 

p^M-S|-l*.^MP (2.26) 
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where "^hi^) are the coefficients of the in-faUing wave-functions, see equation (|2.25| ). 

Equation (|2.26|) shows that once the background metric is known, the spectral 
density is determined by the coefficient that governs the horizon asymptotic of the 
in-faUing fluctuation wave function. 



2.3 The high-frequency limit 

In most of the following we assume there is an asymptotically AdS region r — )■ 
where 

log6(r) ^ + subleading, /(r) ~ 1 + O(r^), r -)■ 0, (2.27) 

and a horizon aX r = where f{rh) = 0, and f'{rh) and 6(r/j) remain finite. The 
black hole temperature is given by: 

4vrT = -/'(r^). (2.28) 

Later we will generalize some of the results to other backgrounds that do not satisfy 
( p.27|) , that will be discussed separately in Section |2.3.3| . 



We make no particular assumptions on the sub leading terms in equation ( 2.27|) . 
In case these sub leading terms actually vanish sufficiently fast as r — )• 0, then the 
metric is asymptotically AdS in the usual sense. 



In we determined the high-frequency behavior of the spectral density from 
equation ( p.26| ) by approximately solving eqs. ( p.21| ) and ( p.22| ) for large u. This 
can be done using an adaptation of the WKB method, as we review in Appendix ^ 
Here, by high-frequency limit, we mean the limit 70; ^ T, A (where A is the scale 
that drives the breaking of conformal symmetry, and can be linked to a deformation 
by a relevant or a marginally relevant operator in the UV). More precisely, the WKB 
calculation in shows that the high-frequency approximation holds for juTs ^ 1. 
Since ^ l/{y/^T) this implies ^/yuJ ^ T. 

On the other hand, even for finite quark mass Mq, we will always restrict 70; ^ 
Mq. The reason is that, in the opposite limit, the description of the heavy quark we 
are using becomes inconsistent. This is due to the fact that for large frequency of 
the fluctuations the accelerations become large, while we are assuming here that the 
configuration is a small perturbation around a steady v = const motion. For large 
70; 3> Mq the backreaction of the fluctuations cannot be neglected. Another, related 
problem is that for 70; ^ Mq, the boundary on-shell action for the fluctuations 
does not reduce to the kinetic action for a quark of mass Mq (as it is the case in 
the opposite limit, Mq ^ 70;) but rather it gives a non- derivative potential term. 
Thus, the holographic pictures of the heavy quark seems to be invahd for modes of 
arbitrarily high frequency. 
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To summarize, the interesting high-frequency regime, for a given temperature, 
quark mass, and quark velocity, is^: 



1 1 T Mn , , 

— <7^<— ^ — — . 2.29 



In this regime, the spectral densities are approximated by 

_2 7^ s r1Juj)Rtp{u)) 

' ^ ^ 27r2£2 1 + (7u;rQ)2 + O ((7u;rQ)4) V ; 

The derivation of this expression is presented in in Appendix |^. Here, rtp{u) is 
the classical turning point of the differential equation for the fluctuations, Rtp{oj) = 
R{rtp{u)), with the function R{r) defined in equation (|2.23|) . This general expression 



simplifies when we take into account the explicit asymptotic form of R{r) and when 
we take Mq — )■ oo, as we will see below. Here, we make a few important remarks. 



The expression ( p.30|) gives the leading high-frequency behavior of the spectral 



densities in terms of the background metric. More specifically, the result is 
completely determined by the UV geometry close to the boundary. The only 
assumptions needed to derive them are that, close to the boundary, the metric 
obeys relaxed AdS asymptotics: 

i rh' 

b(r) -h(r), with ^ as r 0. (2.31) 

r h 

Under this assumption, the classical turning point of the differential equation 
is approximately, 

/2 

Tip ~ — < r, (2.32) 

70; 

for large u. Also, the boundary endpoint rg scales approximately as 1/Mq for 
large^ Mq. Thus for large Mq and large u, the metric functions appearing in 
equation (|2.30| ) are evaluated at a point close to the AdS boundary r = 0. 



For Vfp and tq close to the AdS boundary, one can approximate 



which easily follows from equation ( p.23[ ). 



■*For this interval to be nonempty one needs rg < r^, otherwise the trailing string picture is 
inconsistent even at the level of the background solution. 

^Large compared to the temperature scale, and eventually the scale that breaks conformal in- 
variance in the model. 
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• In the limit Mq — >■ oo, i.e. when the quark is considered non-dynamical, the 
trailing string description is valid for arbitrarily large frequency. We can then 
drop the denominator, and replace (from equations ( p.32| - |2.3^ )). r^pR{rtp) ~ 
rlb\n,)/^^fh\V2hoo), 

Then, equation ( p.3CI| ) simlpifies to: 

p^{u) ^ 7"V||(c^) ^ h'iV2/^u). (2.34) 

This equation shows that, in asymptotically AdS backgrounds, and for Mq ^ •yu, 
the spectral density behaves approximately as at large frequency, as anticipated. 
The leading cubic behavior is deformed by the extra w-dependence induced by 
h{l/uj), which due to the condition (|2.31|) , has a milder dependence than any power- 
law for large u. The function h{u) encodes the deviation from conformal invariance, 
and it may also contain a temperature dependence. 

In the following Subsections we will describe two explicit examples: the conformal 
case (i.e. a pure AdS'-Schwarzschild background), and the logarithmically-deformed 
AdS asymptotics one finds in Einstein-Scalar theories like the Improved Holographic 



QCD setup discussed in This last example captures also the cases that involve 
relevant perturbations in the UV. Finally, we give the generalization of equation 
( |2.34 ) for backgrounds that exhibit scaling but are not asymptotically AdS. 



2.3.1 Conformal diffusion 

When the dual field theory is conformally invariant, the 5D background is an AdS 
black hole, whose metric is given by equation ( |2.11D with: 



6(r) = -, /(r) = 1 - (7^Tr)^ (2.35) 
r 

where £ is the AdS length and T the bulk Hawking temperatures. For a boundary 
quark moving at speed v, the world-sheet black hole temperature T, and the functions 
Z defined in equation (|2.2CI| ) are respectively: 



T, = Z = -. (2.36) 

a/7 7 



Since Z is independent of r, the transverse and longitudinal fluctuations obey the 
same equation (see equations ( ^.21| - P^^2D , and the only difference between p-*- and 
comes from the boundary normalization, i.e. p" = 7^P"'". 

The quark mass Mq determines the boundary endpoint of the trailing string 
through the relation: 
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The fluctuation equation for both transverse and longitudinal modes only de- 
pends on the dimensionless variables x = r /vg and Co = = 



dx 



A 



x'^ ^ ^ , _ , (4a; 



-dj^ix^u] 



X'^ 



+ ^21 _^4^ ^(^>^) = 0- (2-38) 



Let us now specify the general results (|2.30|) to the present situation. In the 
UV limit r — )■ 0, we have R ~ '~f~^i'^ /r"^, and h{r) = 1. From (|2.30|) we obtain the 
following large-cu asymptotics: 

p^iu) = 7-VM ^ (2-39) 
In the equations above, one can substitute the A/" = 4 relation between ig and i, 

The result ( |2.39|) , to leading order in the frequency, will hold in any geometry 
which is asymptotically AdS in the strict sense, i.e. any geometry in which the 



scale factor in the string frame behaves strictly as i/r for r — > 0, i.e. as in ( p. 31 
with h{r) = 1 + 0{r) for r — > (In Einstein-dilaton models, this means that the 
dilaton has to be constant in the UV if the Eintein frame scale factor is asymptotically 
AdS). In these geometries the deviation from conformality will give rise to subleading 
corrections in inverse powers of u with respect to the result ( |2.39| ). In the next 



Subsection we will consider an example where this is not the case, because of the 
presence of a non trivial dilaton (hence h does not go to 1 in the UV), and breaking 
of conformal invariance in the UV deforms the leading behavior. 

2.3.2 Running dilaton backgrounds 

A setup which exhibits breaking of conformal invariance in the UV in a manner that 
mimics QCD can be constructed from the 5D Einstein-dilaton theory with action in 



the Einstein frame p3| , |3^ : 



S = -MlNl I 



4(VA) 



2 



3 A 



2 



(2.40) 



In the holographic interpretation of these models, the scalar A is dual to the running 
coupling Aj of the four- dimensional gauge theory. 

This type of actions provides a general class of holographic dynamics, describing 
a CFT perturbed by a relevant scalar operator. In the conventional case of relevant 
perturbations, setting = log A and adjusting the UV fixed point to be at = 0, 
the near UV potential behaves as V = ^ [1 + C(0^)]. 

For a theory with a UV fixed point at A = 0, the potential should have a regular 
expansion as A — )■ 0, with 

V{X)^^{l + voX + ...) . (2.41) 
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With these requirements, the solutions in the Einstein frame are an asymptotically 
AdS metric, with AdS length i, and a non-trivial profile A(r) which vanishes loga- 
rithmically at the AdS boundary r = 0: 



bEir) 



logr 



9 



8t;ologr \ log r / 

(2.42) 

The scale factor in the above equation is written in the Einstein frame. On the other 



hand, the geometry (|2.11| ) felt by the trailing string is the one in the string frame, 
which in a non-trivial dilaton background is given by 



b{r) 



(2.43) 



Thus the string frame metric is not asymptotically AdS, but it nevertheless obeys 
the relaxed requirement ( p.31[ ), with 

Mr) = X'^'ir) [1 + 0(A(r))] ~ (-^^) [1 + • • •] (2.44) 

For an appropriate choice of the potential ^(A), the models with action ( 2.4C ) and 
asymptotics as in (|2.41|) provide a good holographic dual to large- A^^c 4-dimensional 
pure Yang-Mills theory, at zero and finite temperature pc 



3c 



For a short review of the main features of these models, the reader is referred to 



We now consider the Langevin spectral densities in this class of models. Using 
equation ( p.44|) we obtain from ( p.34|) the large-w behavior: 

01 



Iv^) 27tHI (loga;)4/3 



1 + 



1 



logo; 



(2.45) 



As for the conformal case, we remark that the above results are temperature inde- 
pendent. 

The leading result ~ co^ / (logco)^^^ in (|2.45|) is corrected by a power series in 
Xtp (logcj)"^, which stems from the fact that in deriving equation (|2.45| ) we have 



ignored terms in the fluctuation equation that are proportional to rX'/X and 
also from the series expansion in h{r). From the background Einstein equation these 
terms can be expressed in a power series in A using the superpotential function W{X) 
(see 0): 

oo 



rA7A = AW^(A) = ^PV„A' 



(2.46) 



n=0 



where the coefficients in the power series are determined by the series expansion 
coefficients of V^(A). Thus, the leading result (p.45|) will have the form of a series in 
Xtp (which is a small parameter for large u), i.e. 



log^/3 



1 + E 



^log"a; 



(2.47) 
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The coefficients c„ can in principle be determined order by order from the series 
expansion of the superpotentiaL These coefficients are temperature-independent, as 
shown in Thus, although the leading term is corrected by an infinite power 
series in inverse logarithms, the temperature dependence cannot enter at any order 
in this series, and it is thus confined to smaller power- law corrections in uj. 

2.3.3 Scaling solutions 

The result ([SP can be generalized to solutions that are not asymptotically AdS 
black holes, but are such that the string frame scale factor and blackness function 
behave as simple power laws: 



b{r) 



fir) 



1 



(2.48) 



These backgrounds are solutions of the Einstein-Dilaton system in p + 1 dimensions, 
if the scalar potential is a simple exponential (see Appendix y for details): 



S 



-9 



R 



V{(p) = 2Ae- 



(2.49) 



The interest in these backgrounds is that they arise in the IR behavior of all 
Einstein-dilaton theories with a potential that asymptotes to a single exponential 
at large 0, [^. They also describe non-conformal branes in string theory, and 
their holographic renormalization is well understood |Q. They can be obtained by 
dimensional reduction of AdS black holes in a higher- dimensional pure gravity theory, 
and this explains their scaling properties in terms of a "hidden" higher dimensional 
conformal symmetry . 



The second order equation satisfied by the fluctuation is 



X 



-2a 



v/(l-a:^)(l-(l-t;2) 



2rr4:a 



(2.50) 



v/(l - x^){l - (1 - v^)x^ - v^x^^ 



with 



cil 



V ) 2c U) 

T 



, (2.51) 

It is derived in appendix p. x is a scaled radial coordinate and the various constants 
appearing in the equation are defined in the appendix 0. 

The high-frequency analysis of the world-sheet fluctuations in these backgrounds 
leads to the following spectral density (as we show in detail in Appendix 0): 



7 V" 



e2a 



7 



2a 



7r£2 22^+ir2(a + 1/2) 



2a+l 



[1 + O {{^oorQ)')] . (2.52) 
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The case a = 1 (which arises when 5 = in equation ( |2.49 ) correspond to the AdS 
scahng ~ u^, and (|2.52|) reduces to the conformal resuh (|2.39|) . 



The results in this section all lead to spectral densities that have a divergent 
large-cu behavior, in particular p{u) oc~ for conformal or asymptotically confor- 
mal theories. As discussed in Section 2.1, this cannot be the correct behavior for the 
physical spectral densities that are used to compute real-time Langevin correlators, 
since they cannot satisfy the appropriate dispersion relations ( |2.8| ). For these rela- 
tions to hold, p{u}) should have a well defined Fourier transform, and it should decay 
at least as 1/u at large frequencies. 

To put it differently, the spectral densities we found above would imply that 
dissipative and retardation effects would be stronger for very high energy modes. On 
the other hand, on physical grounds we should expect that the high-frequency modes 
would behave like in the vacuum, i.e. obey a classical equation with no retardation 
and no noise term. 

As we will show in the next Section, the resolution of these problems becomes 
clear once we realize that, even in the zero-temperature vacuum, the spectral den- 
sities computed holographically are non-zero, and display the same large- frequency 
behavior. This suggest that spectral densities satisfying physical requirement at any 
temperature T can be obtained by subtracting the vacuum contribution. In Section 
^ we will show that this is indeed the correct prescription. Before doing that, we 
need to compute the vacuum spectral densities. 



3. The vacuum force correlators 

The zero temperature correlators for an infinitely massive quark moving in the con- 
formal plasma have been analytically derived in |jl^. In this Section we generalize 



those results to a generic T = background with metric given by ( |2.11| ), but with 
/(r) = 1. We will also compute the spectral density in the case of a quark with 
finite (but large) mass, Mq ^ joj (as explained in Section p.3| ), both in the confor- 
mal background and in non conformal backgrounds, such as the scaling backgrounds 
defined by ( |X1^ ). 

In this paper, we restrict to backgrounds such that the string frame scale factor 
b{r) is a monotonically decreasing function of r (from the UV to the IR). This means 
that the gauge theory side at zero-temperature is not confining. In this situation, 
there is no conceptual problem in describing the propagation of a single quark through 
the zero-temperature medium. On the other hand, for backgrounds dual to confining 
gauge theories, this is not evident, and one runs into subtleties. As we will see in 



Section 3.5, these subtleties will result in some immediate problems if we simply try 



to generalize the results obtained in this Section to the confining case. The proper 
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treatment of the dressed Langevin dynamics in the confining case will be given in a 
separate work |]65[| . 

Taking the trailing string anzatz ( p.l3| ), the configuration that minimizes the 
Nambu-Goto action ( p. 12 ) at zero temperature is a straight string stretching from 
the boundary to the interior of the bulk, characterized by ^'{r) = 0. 

Considering the fiuctuating trailing string described by the ansatz ( |2.16| ), we 
can compute the action for the world-sheet fiuctuations. Using the trailing string 
embedding, ^'(r) = 0, the Nambu-Goto action at second order in the fiuctuations 
reads 

2 

gfdjx^dpsx^ + g'fJ2^o.sx''^p6X'' , (3.1; 



S2 = — - / dt dr 



a=l 



where Q^'^ and Q'^^ defined as in ( |2.18|) with Z{r) = I/7 and 

H''f' = ^birfdmg(^-l,^y (3.2) 

Hence the fiuctuation equation at zero temperature becomes 

b' 

^" + 2-^' + 7V* = 0, (3.3) 


both for transverse and longitudinal modes, \E' = e*'^*5X*, i = 1,2, 3. 

In the UV region r ~ this fiuctuation equation has the same form as the one 
for finite temperature. On the other hand, in the IR there is no horizon, and the 
equation depends on the behavior of the scale factor at large r. We will restrict the 
attention to models for which the conformal coordinates extends to infinity, and such 
that the (string frame) scale factor has a power-like behavior in as r — t- +00: 

6(r)~r"'', r -> +00, a>0 (3.4) 

In particular, the conformal case corresponds to a = 1, whereas in the Einstein- 
dilaton models of Section ( p.3.2|) the exponent a is determined by the behavior of 
the dilaton potential for large A. A power-law fall-off of the scale factor as in ( |3.4| ) 
corresponds to a dilaton potential that at large A behaves as V^(A) ~ A'^ with Q < 
4/3. With this parametrization the power a is determined by the relations a = 
(1 - 9gVl6)"^ if log A is not the string dilaton, or a = (1 + SQ/A)-^ if it is [0. 

Using ( p.4|) in the wave equation, we immediately obtain the following asymptotic 
expression for the in-falling wave function 

*«(r,u;)^C^rV'^^ r ^ 00, (3.5) 

As in Section 2, evaluating the conserved flux lm'$*Q''^dr'^ at the horizon we 
obtain the spectral densities: 

^-L _ ^,-2^11 _ k"^(^)]r^oo , ,1^ |2 /q c\ 

Po -1 Po — 2v]-2£2 ^\'~'^\ I'^-t'J 
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The numerator is finite thanks to the asymptotic behavior (|3.4| ). Notice that, 
although the black hole scale factor b{r) reduces to the zero-temperature scale factor 
as — )■ oo, equation ( |3.6| ) is not the limit of equation ( |2.26| ) for — )■ oo, which in 
fact vanishes as ~ r^a. 

We can use the WKB computation outlined in Appendix ^ in order to extract 
the high-frequency behavior of the zero temperature spectral densities. The only 
difference in the WKB computation at zero temperature is in the form of the in- 
falling wave-function, which in this case is given by ( p.5| ). All other steps in the 
determinations of the coefficient \E'o(w) are unchanged, and we arrive at exactly the 
same expression as ( |2.45| ). It now contains zero-temperature background quantities. 
With /(r) = 1, i?(r) = b'^{r)/j, the spectral densities at zero temperature read, in 
the large u regime: 

Notice that the behavior is asymptotically independently of the power-law 
( |3.4|) taken on by the scale factor in the IR. In other words, the high-frequency limit, 
as expected, depends only on geometry in the UV region. Moreover, the leading 
coefficient is exactly the same as at any non-zero T, cfr. equation ( p.30| ). 

Next, we evaluate the explicit large u behavior in the examples of interest. 

3.1 The confer mal case 

In the conformal case, one can solve equation ( |3.3| ) exactly, and evaluate equation 
( ^.24 ) directly in the UV, bypassing the expression ( |3.7| ). The exact solution for ^(r) 



reads: 

\E'(r) = Ai [cos(7Ci;r) + 7Ci;r sin(7Ci;r)] + A2 [sin(7a;r) — 7a;r cos(7a;r)] , (3.8) 
and the spectral densities are given in terms of Ai, A2 by: 

(3.9) 



boundary 



We can immediately deduce the form of the exact spectral densities by setting the 
relation between Ai and A2 from the in-falling condition at infinity ( p.5|) : 



= -iA2. (3.10) 
The normalization at the cutoff boundary r = tq implies 



1 



\M' = \A2? = ,2 - (3.11) 

1 + [l^rq) 
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From (p.9|) we get the exact result: 



7 



2 • 



(3.12) 



In the hmit of infinite mass Mq — )■ oo, or tq — )■ 0, this reproduces the formula for 
the spectral densities in |]T1|, once we use -^^/^^ = V^a^=4. 



3.2 Running dilaton backgrounds 



Speciahsing equation ( p. 7]) to the running dilaton models of Section |2.3.2| , we obtain 
for the zero-temperature spectral densities in the large w, infinite mass limit the same 
expression as in (|2.45| ) 



t'pH^) 



4/3 



27r2 



1 + 



1 



logw 



(3.13) 



(logw)4/3 

3.3 Scaling backrgounds 

For zero-temperature scaling backgrounds, described by equation (|2.48|) with / = 1, 
we can once again solve the vacuum fiuctuation equations ( p. 3D exactly, in terms of 
Bessel functions: 



•qj(r) 



J 



(3.14) 



where Jy and Y^, are the Bessel function of the first and second kind respectively and 
the linear combination is chosen so that the wave function in the IR is an in-falling 
wave. Imposing unit normalization at the boundary r = rq, with joorQ <^ 1, we 
obtain the expression for the spectral density: 



Po 



i-'pI 



fj2a 



^2a^2a+l 



27r2 



|2 -I , hu^rgY ' 
' 2a-l 



(3.15) 



for a> 1/2. Here Ca is the complex coefficient appearing in the leading order in the 
expansion of the Bessel function of order v = —a — 1/2: 



\Ca\ 



22a 
TT 



—TH-a + - 



(3.16) 



The approximation in ( ^.15| ) is valid as long as "^ojrQ ^1. It can be made more 
accurate by adding the next term, O (w^""''^), coming from the Bessel expansion. This 
term becomes more and more important as we send a to 1/2. For half- integer values 
of a a logarithm appears in this term, making it of order {^uirof""^^ logi^uro). 
The exact expression for the spectral density for all values of a > reads 



r'pl 



p2a 



(3.17) 
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For a < some pathologies can arise, as we will describe in Section p.5|. This 



class of models will be further analyzed in a forthcoming paper ||65 



For < a < 1/2, the finite mass correction in the denominator of ( |3.15| ) is slower 
than {■jcorgY. More precisely, it would give a lower power behavior (or power times 
log, for a = 1/2), namely {'jurg)'^"-'^^ . 

The results above, equation ( |3.15D or ( |3.17| ), reduce to the conformal expressions 
of the previous Section if we substitute a = 1. 



As could have been expected, the leading UV asymptotics of the spectral func- 
tions is the same at finite and zero-temperature, and both display a growth as (or 
^"^0-+!^ for scaling solutions). This suggests that we can obtain correlators that are 
well behaved in the UV, simply by subtracting the zero-temperature spectral densi- 
ties from the finite temperature one. This method has been already applied in the 



holographic context for the bulk spectral densities [§5|. In the next section we will 
provide an a priori justification that this is in fact the correct prescription. There 
are two points however that we need to be careful about, if we want to perform such 
a subtraction: 

1. The T = subtraction must kill all terms that grow like and u, and leave at 
most an behavior at high frequency. In order for this to work, all divergent 
terms must have temperature-independent coefficients. 

2. Subtracting the T = correlator should not modify the small frequency behav- 
ior at finite T, which controls the diffusion constants: it would be suspicious 
if the need to regulate a short- distance behavior would also change the long- 
distance physics. 

The first point will be analyzed in detail in Section 4, where we explicitly construct 
the subtracted correlators. The second point is the subject of the next Subsection. 

3.4 The diffusion constant at zero temperature 

In this Subsection we will explicitly show that the zero temperature diffusion constant 
defined by 

.^ImG^ (3.18) 



vanishes for all backgrounds with asymptotics given by ( p^ . As a particular case, 
we recover the vanishing of the zero temperature diffusion constant in the confor- 
mal background, i.e. a = 1. At finite temperature, the computation leading to the 



Langevin diffusion constant was performed in detail in the Appendix B of ||21| , mak- 
ing use of the formula (|2.26|) . Since there is no horizon at T = 0, the method must 
be modified, evaluating the wave function at r — )■ oo, rather than r = r^. 
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First of all, we derive the expression for the diffusion constant in terms of the 
coefficient appearing in the IR behavior of the wave function, ( p.5|) . Substituting 
the solution ( p.5| ) into the definition of the retarded correlator ( |2.24| ), and finally 
using ( p.l8| ), we immediately obtain 



V 



lim \CJ^. 



(3.19) 



The next step is to compute C^. This can be done in analogy with the compu- 
tation of the coefficient \E'/i defined in (|2.25|) , and appearing in ( p.26| ), in the finite 
temperature case (see 1^). Determining requires to match the IR and UV solu- 
tions of ( p.3| ). in the small frequency limit. 

In the IR, b{r) is given by ( |3.4| ) and the solution is written in terms of Bessel 
functions, 



■q)(r) 



" 2 



-e '^^C^^/jLorr" 



C^(caW"" + 4w^+Vi+2") 



u ^0 
r — !■ oo 



(3.20) 



The (UV normalized) zero frequency solution of (|3.3|) on the other hand is given 



by 



^a;=o(r) =Cs + C^ 



dr b 



-2 



Cs + J 



Cv (^l+2a 
+2a V 



+ K) 



r — i- oo 
r 0. 



(3.21) 



Here Ca,da and K are finite w- independent constants. The w-independent constants 
Cs and can be related to if we match the low-frequency behavior of equation 
( p.20|) to the IR behavior of equation ( p.21| ). The outcome is: 



Cs = lim C^ u-'' {ca - daKu'-^'') , C, = (1 + 2a)4 hm C^ 



l+a 



(3.22) 



Imposing the boundary normalization of the zero-frequency wave function, taking its 
asymptotics from equation (|3.21| ), at the cutoff r = rg, implies 



Cs + C^Q = lim C^co " 



Ca ~\~ da 



l + 2a 



K CO 



l+2a 



Hence, at low frequencies, the coefficient C^ evaluates to 



Ca ~\~ da 



1 + 2a 



r?,-K]io 



,l+2a 



~ — , w 0. 

Ca 



(3.23) 



(3.24) 



Substituting ( p.24| ) into the equation for r], ( p.l9| ), brings to the following result 
for the zero temperature diffusion constant: 



ri ~ lim u 



2a 



(3.25) 
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From ( 3.25 ) we conclude that the zero temperature diffusion constant vanishes 
for all backgrounds such that the metric in the IR satisfies (|3^ ), with a > 0. Such 
backgrounds include the conformal metric (a = 1) and all the scaling solutions (a = 
a), for which the result ( p.25| ) agrees with the zero-frequency limit of the exact 
computations in Sections ( |3ri| ) and (p.3|). 

Thus, if we regulate the high-w behavior by subtracting the zero-temperature 
correlator, this does not affect the low- frequency limit. 

Notice in particular that the small frequency limit of the correlator is completely 
determined by the IR part of the geometry, as one might have expected. 

3.5 Extension to confining backgrounds 

We could perform the same analysis that brought to equations ( |3.7| ) and ( p.25| ) 
for backgrounds that, at T = 0, are dual to confining gauge theories. These are 
characterised by a string-frame scale factor that grows in the IR, instead of decreasing 
to zero monotonically. The reason for this is that a necessary and sufficient condition 
for (holographic) confinement is that the scale factor has a minimum at a non-zero 
value at some point in the bulk [^. I In 5D Einstein-dilaton theories with action 
( p. 40 ), confinement requires a potential that grows at least as (see [34| for details): 



y(A) ~ A^/3(logA)^, P>0, A^oo. (3.26) 

Confining backgrounds for which the r-coordinate extends to infinity are character- 
ized < P < 1. For these backgrounds we have a power-law behavior for the string 
frame scale factor in the IR^, like in equation ( p.4| ), of the form: 



6(r)~r-^ a = -j^<0 (3.27) 

The case of the Improved Holographic QCD model corresponds to a = —1/2. 

The analysis of the correlators proceeds like in the previous sections. In partic- 
ular, the high-frequency limit is still given by the universal formula ( p.7| ), since it is 
only sensitive to the UV. However, the analysis of the low-frequency limit for a < 
gives a result that at first sight seems puzzling: 

77 ~ lim (3.28) 

with a log^cu correction in the special case a = —1/2. 

The problem with this behavior is that, for the range — 1 < a < (which in fact 
contains the Improved Holographic QCD model), the diffusion constant diverges at 
T = 0. This suggests that the appropriate string solution is problematic in the IR. 

A resolution of this puzzle lies in the fact that in such cases another embedding 
solution is relevant. This will be analyzed in a separate work [|65 |. 



^Note that the Einstein frame scale factor behaves exponentiahy in the same regime logbE 

^2a+l 
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4. The dressed Langevin evolution 



Having analyzed the finite- and zero-temperature Langevin correlators in the previous 
sections, we proceed to show that a natural choice for a UV-regulated propagator 
Gr{uj) that has the correct fall-off at large u is given, in the limit of Mq/u 3> 1, by 
simply defining: 

Gn{co) = Gnico)-Go{co). (4.1) 

From the previous Sections we know that this procedure is consistent for all 
non-confining and scaling backgrounds, since it evidently takes care of the leading 
divergent coefficient, and it does not modify the low-frequency behavior of the finite- 
T correlators. At the end of this section we will also show that, after the subtraction, 
the dressed correlator has the correct fall-off at large frequency to be consistent with 
a regular short-time evolution. 

Here, by "zero-temperature" correlator, we mean a correlator evaluated in the 
background whose metric coincides with the vacuum metric, but it is heated up at 
a temperature T, i.e. the thermal gas background. This is necessary if we want to 
stay in the framework of the canonical ensemble, where T is a fixed parameter. The 
thermal gas is a solution of Einstein equation that corresponds to a saddle point of 
the action, the true minimum being the black hole of the same temperature. 

4.1 Path integral treatment 

A physical correlated Langevin diffusion has force correlators that are regular func- 
tions of the time difference. In particular, as we have seen in Section ^, an application 
of the formalism at T = (in AdS) would provide a classical Langevin evolution (|2.2| ) 
with a friction term analogous to the retarded correlator computed in AdS, that as 
we have seen is non-trivial and UV divergent. However, we do not expect that a 
(free) heavy quark feels a force in the vacuum. This suggests that the proper phys- 
ical coordinates of the heavy quark should be defined with respect to the T = 
background. We will do it explicitly below, using the path integral formalism. 

We start from the Schwinger-Keldysh path integral, following the treatment 
found in ITB 

VXlVXr e^^«(^«)-^^^(^^), (4.2) 



with X* = Xq + 6X^ where Xq describes the classical motion of the dragging string 
and 6X^ is the fluctuation and 

Six) = SiXo) + J d^^P^SX' -\j d'^i GtfiXo)dJX%5X^ + OUSXf). (4.3) 

In this language, the trailing string equations are 

do^Pt = (4.4) 
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and the classical equations for the fluctuations are 



d^[Gtf{Xo)dp5Xn=0. 



It is convenient to Fourier transform 



du 



/vr 



(4.5) 



(4.6) 



We can split the fields 5Xi into a boundary piece Xi and the bulk piece 5X^. Inte- 
grating out the bulk pieces in ( ^4.2|) we obtain the boundary action as 



(4.7) 



with 



iSh = —i 



^C(-^)Gg(^)C(u;)-^ / ^Ca{-u^)G%^{^)a{co) 
du 



2tt 



(4.8) 



where we passed to the advanced and retarded basis, and Gr, Gsym are the standard 
retarded and symmetric correlators, and C*('^) = Ci^^f^b) are the boundary values 
of the Fourier space variables. 
We now use the identity 



exp 



du 



2 ./ 27r 

r \ 1 f du 

X J VC exp j - 

The path integral becomes 



a'^)G%„,{u)ci{u 



G-L\ 



det ( j X 



(4.9) 



du 
2^' 



dfj, e" 



with 



and 



(4.10) 



(4.11) 



S 



du 
2^ 



(4.12) 

Integrating out Q gives a functional 5-function that imposes the Langevin equa- 
tion, which can be written in configuration space as 



dt'G%{t,t')ci{t') + p\t) = e{t). 



(4.13) 
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The noise variable has two-point function G^^y^. 

The previous derivation holds also in the case of a quark propagating in the 
vacuum, since even in pure AdS at zero temperature (and more generally in the 
non-conformal vacuum solutions) Ggym and Gr are non-trivial. In particular, in the 
vacuum solution we would get a Langevin equation like ( [4.13|) with Gr replaced by 



G^, the vacuum retarded correlator, and with noise distributed according go G°y^. 

The fact that a particle in vacuum seems to undergo dissipation induced by 
quantum vacuum fluctuations seems rather unphysical. This is an indication that 
we need to modify the definition of the path integral in order to make sure that 
the average coordinate Cr('^) of a isolated particle in vacuum evolves according to 
its classical equation of motion, P{C,r{^)) = 0. For simplicity we will assume that 
there is no external potential, thus the classical equation of motion is a linear term 
in the average coordinate: P* = JH^^ {co)(rj{^^) (e-g- for a non-relativistic particle 
M'^iu) = Mqu^5'^ ). 

We now proceed to define a path integral for a "dressed" string variable in such 
a way as to obtain the classical trajectory in the vacuum. We give this redefinition 
at the level of the boundary variables, i.e. after the bulk part of the string has been 
integrated out. The change in the path integral is composed of two parts: 

1. A change of variables: 

Ciu^)=[iM{co) + GUco)r'M{co)]lC^ico), C(c^) = CM- (4.14) 

We will take the dressed variables (r and (a as the true physical variables that 
describe the boundary quark. 

2. A modified integration measure. We redefine the integration measure of the 
path integral to be used in ( [4.7D by: 



dfl = VCa'DCr exp 



1 1 



(4.15) 



where VC^ = Ylu^^^^i^) bare functional measure. 

We define the dressed boundary path integral at zero temperature as: 

Zo = j dfiexp'^°^^^''^^^ (4.16) 

Using the explicit form of the action from ( [4.8| ) (with the zero-temperature Gr and 
Gsym) we get immediately: 

Z, = j VCM^e^f mM%: jy^^s (^Mt cm) , (4.17) 
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where A/" is a normahzation constant. In other words, the vacuum path integral is 
concentrated on the classical trajectory with no dissipation for (r, as required. This 
justifies our interpretation of the "dressed" variables ( as the physical ones describing 
the motion of the heavy quark. 

Now we come to the non-trivial part, i.e. we define in the same way the path 
integral at finite temperature: 



dfi exp 



(4.18) 



Performing the change of variables ( ^.141 ) and using the modified integration measure 
( [4. 151 ), it becomes: 



Z 



[m+Gr)[m+G''j^ 



M 



(4.19) 



Following the same procedure that led to equation ( [4.13|) , we obtain a Langevin 
dynamics for the dressed variable C,a defined by the dressed Langevin correlators: 



{Gr)\= \[M + Gn)[M + G"RY'M 



Mi 



G 



sym 



ym 



^ sym } j ■ 

(4.20) 

For the symmetric correlator, the dressing simply amounts to the subtraction of 
the zero-temperature term. To obtain a simpler expression for the dressed Gr, we 
first consider the non-relativistic case where A^*(a;) = Mqu'^Sj. In this case, the 
first equation in ( [4.20| ) can be written as: 



{GnV = \\{1 + Gr/Mqu') (1 + GyMguy 



5] \ Mou\ 



(4.21) 



Throughout the paper we have been making the assumption that uj/Mq ^ 1, in 
order for the holographic picture to be reliable. Since the zero temperature correlator 
behaves as both at large and small frequencies, we can expand the denominator 
and arrive at: 

'GrX, = {Gr - G?,);. [1 + {00 /Mq)] . (4.22) 



Thus, the dressed retarded correlator also reduces, in this limit, to the subtracted 
correlator. 

In the relativistic case, some extra work is required, since the effective mass 
matrix is not proportional to the identity, but is given by equation ( p.4|) . Introducing 
the longitudinal and transverse projectors. 



n 



6i 



VjVj 

^,2 ' 



n" 



ViVj 



we can write: 



eff 



(4.23) 
(4.24) 
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Correspondingly, we split the retarded correlator in a longitudinal and transverse 
part, 

Gn = (n^ + nil) GR = Gj, + gI (4.25) 



and similarly for G^. Next, we rewrite the exrpression for Gfi in ( [4.2(J| ) as: 

Gr = {1 + M'^Gr){1 + M-^G"ji)-^M - M (4.26) 
and notice that A^^^G^ is small for small jcjMq: 

Thus we can approximate (1 + Ai^^G°f,)^^ ^ (1 — M.^^G'}^), and after some algebra 
equation ( [4. 261 ) becomes: 

Gr ~ (Gr -G'i,)M = [Gr - G° )^ + (Gr - G^)" = (Gr - G%). (4.28) 

Thus, with the assumption u ^ Mq we have the desired result that the dressed 
correlators are simply obtained by subtraction of the zero temperature correlator both 
for the retarded and the symetric Green's functions: 

Gr — Gr — G^, Ggym = Gsym ~ G°ym- (4.29) 

The correlators G r and Ggym are the ones that should enter the physical Langevin 
diffusion in the plasma. Gr defines the physical spectral densities, 

p{uj) = p{u)~p°{u) (4.30) 

To complete the discussion, we establish the relation between Gr{uj) and Gsym{u}). 
For the bare correlator in the black hole background we have the Einstein relations 
at temperature T = Tg. It is reasonable to chose the temperature of the thermal gas 
background as the same Ts, such that two trailing string live in the same ensemble. 
Thus we have: 



- s 



GsyU^) = coth ( ^ ) Im Gr{uj), G°yjuj) 



coth (^^^ Im G?j(a;) (4.31) 



If we make this choice of ensemble, it follows that the dressed correlators satisfies 
the standard Einstein relation: 

Gsym = coth (^^^ Im Gr. (4.32) 
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4.2 UV behavior of the dressed spectral densities 

Having defined the physical Langevin spectral densities p{u) at finite T by subtract- 
ing the zero-temperature counterparts, our final task is to show that they display 
the correct large-w fall off, at least as 

As it is evident from Sections 2 and 3, p{u) does not contain the behavior at 
large frequency, since it cancels in the subtraction. However, this is not enough, as it 
could be that some terms that diverge linearly in u are left over after the subtraction, 
which would be unacceptable. This is not the case, as we show below. 

Let us recall that the regime we are interested in is the high-frequency limit, 
70; ^ 1/ts, and at the same time we need jojrQ ^ 1 to be in the regime where the 
quark can be considered non-dynamical. 

First, let us analyze the temperature-dependence of (|2.30|) . For this, we need the 



subleading behavior of the background functions entering the fiuctuation equations 
near the boundary. The crucial point is that temperature dependence manifests 
itself, in asymptotically AdS backgrounds, as terms which all vanish as r"^ close to 
the boundary. 

For asymptotically AdS black holes, as r — )■ 0, we have 

/(r)~l-C/ (4.33) 

where C is a constant proportional to the product ST of the black hole. Thus, from 
the expression (|2.23| ) for R{r), we can write: 



R(r) 



(4.34) 



7 

We also know, by a slight generalization of the discussion in [P7|, that close to the 



boundary'': 

b{r) = boir)[l + g{T)g,{ry], (4.35) 

where g{T) is a temperature-dependent constant which controls the conformal anomaly, 
and gi{r) is a slowly- varying function that does not modify the power- law (other- 
wise consistency of Einstein's equations would spoil the AdS black hole asymptotics), 
but whose exact form depends on the specific model, and on whether the Einstein 
and string frame scale factor coincide. 

The effective Schrodinger potential for the wave-function is given in equation 
( |A.4D . From the discussion above, it is clear that the only temperature-dependent 



corrections to the leading behavior, close to the boundary, scale as: 

SV{r) ~ r\ (4.36) 



''In general the scale factor will be temperature dependent. The only case when this is not so is 
that of ^dS'-Schwarzschild black holes 
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This would give rise to corrections to the wave functions that close to the boundary 
vanish as r'^, and when evaluated at tq it can only give rise to a temperature- 
dependence in the UV WKB coefficients which is of order with respect to the 
leading behavior. 

Close to the horizon, on the other hand, the expression ( p.25| ) does not get modi- 
fied by finite l/u corrections. Thus, if we follow the matching of the UV and horizon 
coefficients performed in Appendix |^, we conclude tht the only temperature depen- 
dence in ( |2.30|) enters at subleading 1/(70;)^, and there are no dangerous 1/(70;)^ 
terms. 

We can then write, from the general expressions (|2.30|) and ( p.7|) at large u (for 
e.g. the ± spectral density): 



P (w) 



27r2£2 1 + (^ujrgY 



(4.37) 

where rtp^o and ho{r) are respectively the WKB turning point and the scale factor at 
zero-temperature, the last brackets include the temperature-dependent corrections to 
( p.30|) discussed above, and we have dropped the {'^ojxqY terms in the denominator^. 



The term in the ffist bracket also gives rise to terms that grow no faster than 
l/uj^: for large frequency the turning point is close to the boundary, rtp ~ 1/(70;) ^ 
rg, both at zero and at finite temperature; the functional dependence on rtp scales 
as r^p due to equations ( [4.33[ - P755D , and the turning point is determined by the 
Schrodinger potential V{r), which is also modified at (9(r^). Thus the turning points 
at a fixed frequency 70; are shifted at most by the same fractional amount, 

rtp = rtp,o [1 + 92irtp, T)rfp] , (4.38) 

where again g2{r, T) is another slowly varying function which can be determined from 
the Schrodinger potential. 

Collecting results ( |4.34| - |4.38|) in the expression (|4.37|) we arrive at: 



^ ~ 27tX 1 + il^rgy J^U ' ^ ^ ^ 

where g{u]T,rQ) is a temperature-dependent function, slowly-varying (with respect 
to a power law) in u, which collects all the various contributions discussed above. 
This is the final result, that shows that the dressed correlator we have defined has a 
well-defined Fourier transform, thus it satisfies all the appropriate dispersion relations 
and provides a consistent Langevin process at short times. 



^these terms can easily be evaluated explicitly from the WKB calculation in Appendix For 
example, the full expression in the case of IHQCD can be found in They depend on h(r) 

evaluated at the turning point, and by the argument presented below they cannot contribute more 
than l/o; in the UV to the subtracted spectral densities. 
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5. Examples of dressed Langevin correlators 



In this Section we present the numerical evaluation of the full dressed correlators, as 
a function of u, in two interesting models: the conformal case, and the non-confining 
scaling backgrounds. 

5.1 Dressed correlators in the conformal case 

In the conformal case all temperatures above zero are equivalent. However, since 
we introduce a cutoff determined by the quark mass as in (|2.37|) , the mass-to- 
temperature ratio is related to the ratio rh/rq = y/jrg/rQ (and by the velocity, 
through 7) as follows: 

Mq ^ v/AAr=4 rh , 
T 2 TQ- ^ ■ ^ 

In Figure |l] we show the imaginary part of the retarded correlator, rescaled by 
y/jT^, and as a function of the dimension- less frequency u = 'juvg, for a specific value 
of Mq/T and different values of p/Mg = — 1. On the left, the renormalized 
correlator is compared to the finite temperature and zero temperature correlators 
for fixed momentum, while on the right we show the renormalized correlators for 
different momenta. 

We chose to use the dimension- less frequency u, because it encodes the whole 
dependence of the rescaled correlator on the temperature and on the velocity, in 




(b) 

Figure 1: The imaginary part of the retarded correlator, rescaled by ^y^T"^, is shown 
as a function of the dimension- less frequency jujrs- For both figures we set Vfi/rq = 
17, implying Mq/T = 20. In (a) we plot the numeric result for the finite temperature 
correlator (red plain curve) , the zero temperature analytic result (black dashed curve) and 
the renormalized correlator magnified by a factor 100 (blue dotted curve), for rq/rg = 
0.19, meaning p/Mq = 10. In (b) we consider three different values for the ratio p/Mg, 
p/Mq = 1, 10, 100 (rq/rg = 0.07, 0.19, 0.59), and show the renormalized correlator for each 
of them. 



a 
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the limit of infinitely massive quarks (when the cutoff is sent to zero) . This follows 
from the fact that the fluctuation equations only depends on the frequency, temper- 
ature, velocity and radial coordinate through this combination and through the ratio 
dimension-less coordinate x = r/r^ (as described in Subsection 2.3.1| ). 



However, in the case of a finite cutoff, there appears an additional dependence 
on the velocity, through the ratio tq/ts = xq, leading to different curves for different 
momenta. This stems from the fact that the flux must now be evaluated at r = rg, 
rather that at r = as in the infinite mass case, where the dependence on the velocity 
and temperature is all taken into account by the rescaled correlator and frequency. 
Hence, as we vary the velocity in the finite mass case, we change the value of the 
momentum-to-mass ratio p/Mq, which is related to the ratio tq/ts by 



p 16 (MQXfrQ 



(5.2) 



Moreover, we note that the cutoff rg should be kept fixed if we want to keep the 
quark mass fixed. 

In Figure |l|(a) we note that the imaginary part of the retarded correlator displays 
the expected behavior at high frequencies: cubic at intermediate frequencies (namely 
for frequencies still much lower than l/'yurq) and linear at higher frequencies. The 
renormalized correlator quickly goes to zero as the frequency grows. 

The dressed correlators plotted in Figure ^(b) show the dependence on the ve- 
locity that arises due to the finite mass. The rescaled dressed correlators depend 
both on the velocity and on the temperature only through the ratio rg/r^. 

5.2 Dressed correlators in scaling backgrounds 

The scaling backgrounds described by the metric functions in ( |2.48| ), with a > 0, 
share a similar behavior with the conformal AdS^ model. In this Subsection we 
show the numerical results for the dressed correlator in this class of models. The 
zero temperature correlator is derived analytically, using the exact expression (|3.17| ) , 
while the finite temperature correlator is computed numerically. 

The presence of a finite mass introduces a scale in these models. The mass-to- 
temperature ratio is given in terms of the quantity Vh/rq by 



(5.3) 



(£T)2«-1 C £2 y^^^ 

We restrict to the region of parameters of the scaling backgrounds where d > 1/2, 



in order to produce a consistent spectral density, as discussed in Section The 
plots shown in figure Q represent the imaginary part of the retarded correlator divided 
by a temperature and velocity dependent factor, ^2(2a-i)/c2n2a+i^ which encodes the 
whole dependence on T and v in the infinite mass limit. As in the conformal case. 
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(a) (b) 

Figure 2: The imaginary part of the retarded correlator, rescaled by 'j/2{2a.-i)/c2^2a-i-i^ -g 
shown as a function of the dimension- less frequency jujrg- For both figures we chose the 
background parameters such that a = 5/6 and c = 23/5 inp + 1 = 5 dimensions (this 
means = and = 1/9). We also set r^/rg = 30, implying IMq/iiTf^-^ = 2. 
As for the conformal case, in (a) we plot the numeric result for the finite temperature 
correlator (red plain curve), the zero temperature analytic result (black dashed curve) and 
the renormalized correlator magnified by a factor 100 (blue dotted curve), for rq/rg = 
0.09, meaning p/Mq = 10. In (b) we consider three different values for the ratio p/Mq, 
p/Mq = 1, 10, 100 [rq/rg = 0.04, 0.09, 0.25), and show the renormalized correlator for each 
of them. 



there is an extra dependence on v (at fixed T and Mq) when the mass is finite. This 
dependence arises through the ratio rq/rs, appearing in the O {{jurq)"^) corrections 
- and higher orders - in ( p. 15 ) and ( 3.17 ). 

The momentum-to-mass ratio can be expressed in terms of rq / as 



P 

Mr 



2ti (ii eur 



^rq_ 



- 1. 



2a - 1 £2 (£T)2«-i y \47rrsJ ^' ^^'^'^ 

So at fixed iMq/{iTy'^'^ (or, equivalently, at fixed Vh/rq) the spectral densities 
corresponding to different momenta p/Mq (or different rq/rg) will have a different 
behavior. This is shown explicitly in figure @. 
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Appendix 



A. The high-frequency expansion 



In this appendix we review and expand the result obtained in |^ about the high- 
frequency behavior of the spectral densities. We will perform the detailed analysis 
only for transverse modes. The same method can be obtained for longitudinal modes 
with minor modifications. The analysis we provide here is general, and covers both 
asymptotically AdS solutions, and scaling backgrounds. 
The starting point is the wave equations (|2.21|) : 



a^(i?9,^(r,w)) + ^^(r,u;) = (A.l) 

where 

R{r) = V(/(r)-t;2)(64(r)/(r)-C2), C = vh\r,). (A.2) 
This equation can be put in Schroedinger form with the redefinition: 



= R'/\r)<i!, (A.3) 



in terms of which it reads: 



_^"(^) + V{r)^P{r) = 0, V{r) = + ^(logi?)" + l(logi?)''. (A.4) 

The first term in the potential is the only important one for large u, except in a 
small region close to the boundary, where R ~ 6^/7: since b{r) diverges as a power- 
law for small r, close to the boundary the u independent terms grow approximately 
as r^, whereas the first term stays constant. Thus, for large u, we can neglect the 
w- independent terms everywhere except in the region < r < 1/(70;). Elsewhere, 
one can use the WKB approximation to compute the solution {WKB region). On 
the other hand, it is possible to approximately solve equation ( |A.4| ) analytically for 
< r ^ r^, {UV region), which for large frequency 70; ^ l/r^ overlaps with the 
WKB region. This allows to determine the coefficient of equation ( p.25| ), thus the 
spectral density ( p.26| ) . This procedure is explained in detail in the next Subsection. 



A.l The WKB wavefunctions 

In the range < r < we can identify three regions: 

• WKB region: 1 / {'ju) ^ r < 

• UV region: rq < r <^ rg. 

• Overlap region: 1/(70;) ^ r ^ r^. 
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The overlap region exists for large frequencies, ^ l/r^. Notice that the WKB 
region contains the horizon. The UV boundary r = rg is finite for finite quark mass, 
whereas rg ^ for strictly non-dynamical quarks. 

Below we solve for \l/(r, w) separately in the UV and WKB regions, and match 
the coefficients in the overlap region. 

WKB region 

In this region, we can approximate 

V{r) ^ ^ (A.5) 



In the WKB approximation, the independent solutions of equation ( [A.4|) are 

^ = 1/^1/4 gxp (±1 [ ^/V{r) 



Using the explicit form ( |A.5| ) for the potential, we can write \& = i? ^/"^ip as: 



b(r ) r b'^(r') 

^WKB = ^^-^ J dr'^ (A.6) 

where the normalization has been chosen for later convenience. The choice of a purely 

ingoing wave is the one that gives the correct behavior at the horizon: as r — )■ r^, 
R{r) behaves as: 

R{r) ~ AtcTs 6^(rJ (r^ - r), r -> (A. 7) 
thus we find, close to the horizon, the appropriate in-falling solution: 

^^yi^B^^h(r.-r)-*'^/^'^^^ (A.8) 

UV region 

For small r ^ the effect of temperature is negligible, as / ~ 1. In this region we 
can approximate 

R ~ (A.9) 

and equation ( A.4| ) becomes simply: 



+ y^tp = u^-f^ij. (A.IO) 
V R 

To find approximate solutions of this equation, we distinguish two situations: 

1. power-like case: b{r) ~ {i/ry 

In case b{r) is a pure power law at small r, b{r) ^ {i/ry (as is the case in 
asymptotically AdS backgrounds, or in the scaling backgrounds analysed in 
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Section |2.3.3| , then {y/R)"/y/R ^ p{p + l)/r^ and the solution for ^ = y/R^p is 
simply given in terms of Bessel functions: 

^^^^''^ = -7^^ [Jp+idiujr) + tNp+y2h(^r)] (A.ll) 

This linear combination chosen in such a way that it connects with the ingoing 
wave (|A.6| ) in the overlap region: for r/r^ <^ 1 but •jur 1 we can use the 
large- argument asymptotic of the Bessel functions: 

Ju{x) ~ \/ — sin (x — UTi /2 — 7r/4) , Ny{x) ~ \ — cos (x — vn /2 — 7r/4) 

V TTX V VTX 

(A.12) 



Using this approximation as well as equation ( |A.9|) to compare the two forms 



(|A.11| ) and (|A.6| ) in the overlap region, we find the relation between the coef- 



fifients A and '^h (up to an oveall constant phase): 

A\— = ^,6(r,) (A.13) 

On the other hand A is fixed by requiring unit normalization '^{tq) = 1. For 
small 7wrQ we can use the small argument expansion of the Bessel functions: 

to obtain the normalization (using R ^ — (£/r)^^/7): 



From equations ( A.13 - A.15 ) we determine '^h-, which inserted in equation ( p.26| ) 



determines the spectral density: 

The conformal case of AdS asymptotics corresponds to p = 1. For this value 



of p we recover the results found in [21]. In this case the Bessel functions have 



a simple expression in terms of sines and cosines, and it is not hard to specify 
the first correction to the leading result: 
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2. corrected power-law b{r) 

Let us now take p = 1, but consider the more general case of the asymptotics 
( p.31| ), where h{r) does not approach a constant as r — )■ (otherwise we would 
fall back in the previous case with p = 1). In the following we will restrict 
to approximate AdS asymptotics with p = 1, but the method can be easily 
generalized for any p. 

With the asymptotics ( |2.31| ), the Schroedinger potential in the UV becomes: 

2 / 2rh' yh"\ . , , 
V^~— 1 ^ + ^-r -7V (A.18) 



\ h h ^ 

If the /i-dependent terms in the parenthesys are small for small r, we may think 
that we can still approximate the solution as a linear combination of the form: 

^(r) = [j3/2(7a;r) + iNs/2{lcor)] (A. 19) 

Expanding this for r — )■ leads to ip c:^ Cg/r + C^r"^ for some constants Cg and 
a, which would lead to ^ ~ R-^/\r){CJr + Cy) ~ h-^/^{r){C^ + Cy). 
This is not, however, the correct behavior for small r, which can be obtained 
directly by integrating equation ( |A.1|) for r ^ 1/ (70;) neglecting the cu^ term: 

r dr' 

In order to get a consistent solution in the whole UV region, we must correct 
equation ( |A.19| ) as follows 



^Pir) ~ [A2F2ir,u)Js/2ilujr)+tAiF,{r,u)Ny2ilcur)] (A.21) 

where Fi and F2 are functions that are slowly varying, and behave as: 

Fi ~ yR{r)y/^, F2 ~ yR{r))-^/\ r ^ 0. (A.22) 

With these ansatz, the wavefunction = R'^^^ip has the correct behavior at 
small r, as one can see by expanding the Bessel functions and comparing with 
equation ( |A.20| ). 

On the other hand, as the solution enters the overlap region (r > 1/(70;)) , the 
functions Fi and F2 must approach two constants, in order to match the WKB 
solution. Since they are slowly varying, we can estimate these constants to be 
the values 



Fi(r) -^Fi ~ (r"i?)^/2 ^ F2 ^ F2 - [r'R)-'^'' , r > 1/(70;) 



-1/2 



rtp 



rtp 

(A.23) 
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where rtp is the classical turning point of the Scrhoedinger equation, and it 
marks the "boundary" between the UV region and the WKB region. For large 
a; it is approximately r^p ~ a/2/ (70;) as can be seen from equation ( |A.18| ) Now 
the matching in the overlap region leads to: 



AiFi = A2F2 = <ilhb{r. 



ITU 




(A.24) 



and the unit normalization Psi{rQ) = 1 is: 



1 =iAi 



l-i- 



Ao ( F.J, 



2 -'3/2 



Ax \F1Ny2 



rQ 



(A.25) 



To lowest order in rq, using the asymptotics ( |A.22[ ) and those of the Bessel 
functions, this gives: 



A, = y^i^oof' [1 + (7u;r)2 ^ ^ {(^oorQ)^)] ' . 



(A.26) 



From equation ( [A.23| ) and ( [A.24| ) we can determine which inserted in equa- 
tion ( |2.26| ) gives: 



/ ^ ^7'oo'rlRinp) [1 + i^oorQ)' + O {{loorQ)')]-' 



(A.27) 



A similar procedure gives p"(c<;), the only difference being an extra factor 7^ in 
the numerator. 



A. 2 WKB at zero temperature 

At zero temperature the blackness function /(r) = 1 and there is no horizon. Cor- 
respondingly, we have 

R(r) = (A.28) 

7 

We consider backgrounds in which the conformal coordinate reaches +00 in the 
interior, where the scale factor behaves as a power-law, 

b{r) ~ r"", a > 0, r ^ 00 (A.29) 

The only change with respect to the previous Section is that the WKB region now 
covers the deep IR, and the WKB wavefunction ( [A.6| ) is given by 

"^WKB = C^r'^e'^^ (A.30) 

One then repeats exactly the same steps as in the previous Section in the UV, to 
determine C^, and the result is the same as (|A.27| ), with R{r) replaced by b'^{r)/'j. 
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B. The boundary effective action 



In this appendix we compute the boundary effective action for the world-sheet fluc- 
tuations, at the boundary point tq. We concentrate on the real part of the quadratic 
on shell action, which should reproduce the corresponding action for the fluctuations 
of a relativistic particle moving at constant speed in four dimensions. We will see 
that this is the case only in the limit 'Jootq — )■ 0. 

We restrict our attention to (possibly corrected) AdS asymptotics. Also, we 
only consider here one transverse fluctuation; generalizaztion to longitudinal ones 
is straightforward. We start with the quadratic action for world-sheet fluctuations, 
equation ( p.l7|) . 



1 r+QO pTs 

S^'^ = J drjdr [H^^d^^dp^] . (B.l) 

Evaluating the action on shell should produce boundary terms that, in the limit 
Tq — )■ 0, should match the form of boundary diff-invariant counterterms for the 
source term in the fluctuation, i.e. \Ef(rQ,t). The only such term, in a homogeneous 
medium, is the one coming from kinetic action for a relativistic particle, 

Sbdr = Mq j drsjx^X-h^,, (B.2) 



where hp,v is the induced boundary metric. Expanding ( [B.2|) around a constant speed 
trajectory with a transverse fluctuation, X° = t,X = ^Q\vt\'^{€)u^ with u^-v = 0, 
we obtain (going to Fourier space): 

C,(2) Mq f^°° du 2 2|,T,^ M2 /T3 Q^ 

^... = ^y_ ^7-|vI^(c.)| (B.3) 

Now we evaluate equation the bulk action, ( [B.l| ), on the in-falling solution 
\E'(r, cu), but relaxing the unit boundary normalization, i.e. we take \E'(rQ,ci;) = \E'(a;) 
to be arbitrary. After an integration by parts and use of the field equations for 
\E'(r, w), as well as H^^ = R{r), we obtain: 

S^'^ = 4^ (/ ^^*{co)R{rQ)dr^{co) - J ^^*{r,,co)R{r,)dMrs,u:)^ (B.4) 

The second term arises from the world-sheet horizon, and can be regarded as an IR 
modification to the quark mass. We are interested here in the boundary contribution, 
so we will focus on the first term only. 

In order to evaluate the first term, we use the high-frequency solutions found 
in Appendix |^. We need to separate the fluctuations in two sets, according to the 
value of u. 
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1. u; > l/iiTQ) 

In this case rq belongs to the WKB region, and we can use equation (|A.6|) , 
where now ^>hb{rs) — >■ \l/(a;). Inserting this expression in the first term in 
equation ([B.4| ) we find: 



(a;>l/(7rQ)) 
UV 



duj 
2^ 



b' R 



+ iu 



(B.5) 



Comparing with equation ( p.3|) we see that the real part of this expression does 
not correspond to a boundary kinetic term for the string endpoint as it is of 
order zero in u. Rather, it gives rise to a harmonic potential. 

2. u; < 1/(7^^0) 

In this case, rg is in the UV region, but not in the boundary region. We can 
use \E'(r) = R~^^'^ip{r), with ■ip{r) given by ( |A.21| ), where now the coefficients 
Ai, A2 are arbitrary functions of u. It is convenient to rewrite the wavefunction 
in the form: 



r/rg Fi{r) A^3/2(7C^r) 
R{r)/R{rQ)F,{rQ)Ny2{lujrQ) 



F2{lo) F2{r) Ji,/2{l^r) 
Fi{uj) F^{r) N^/2{lujr) 



(B.6) 

Since we are in the region rq <^ 1/(70;), we can use the approximation ( |A.22| ) 
for the functions Fi and F2: 



r/rq Fi{r) 
R{r)/R{rq)F,{rq) 



[r/rq) 



3/2 



(B.7) 



Moreover, we can use the small argument expansion A'3/2(x) ^ —y^2/'Kx^^^'^ — 
x^/^/-\/27r + . . . both in the numerator and denominator (it is necessary to go 
to first subleading order to obtain a non-zero result): 



^3/2(7^^) 
N3/2{lujrq) 



(rq/r 



1 + ^ + O {^urY 



(B.8) 



The term inside the parenthesys of equation (p.6| ) , on the other hand, is 
approximated by: 



1 - i 



r'^R{r) 
r%R{rtp 



^3^3j,3 



+ O {inujrf) 



(B.9) 



as r — )■ 0, as it results from the expasion ( |A.14| ) and from eqs. ( |A.22f]Xl23| ). 
Putting everything together we obtain, in the limit of small jurq: 



1 + ^ + O {^urf 



[I + lO [inurf)] 



(B.IO) 
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The real part of the UV boundary action in ( p.4|) is therefore: 



(uKl/i-yrQ)) 
UV 



I'^I<1/(7'-q) 



(B.ll) 



This has the same form of a mass term, with the effective quark mass given 
by: 



Mr 



(B.12) 



Notice that this diverges as we remove the cutoff, rg — )■ 0: as it should be, the 
divergence is cancelled precisely by the boundary counterterm ( B.2 ), as was 
noted in Also, notice that in the conformal case we recover the relation 
Mq = (£V27r£2)l/rQ. 



C. The neutral scaling solutions 

We start with the action 



with 



In the domain-wall coordinate system, 

ds2 = [-f{r)dt^ + dx^dx'] + 
we have the scaling solutions 



dr^ 



(C.l) 



(C.2) 



(C.3) 



2p 



/=l-(^)^" , 



The metric can be rewritten by a charge of coordinates 

t 



w = r (p-i)"*^ , t — )■ 



X 



(p-l)52 



X 



1 



(p-l)<52 



as conformal to an AdS-like black hole, 



^ - fde + dx'dx' 
ds^ = e^x _l 



(C.4) 



(C.5) 



(C.6) 
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with 



^■2 



— d 1 I in \ 2-(p-i)5^ 



Prom this we can rewrite the metric in the conformal (and Einstein) frame 



ds' = b{w) 
with 



-fdr + dx'dx' + 



/ 



(C.8) 



For 5 = wc obtain the AdSp^i solution. Note that always 2p — {p — 1)6^ > 
in order to satisfy the Gubser bound. Therefore, for 6 > 0, the sign of all a, c, d 
is always the same. Also, in the regime 5^ < c e [p, oo), while in the regime 

<^'>^,ce(-oo,o] 

We now consider a dragging string moving in this metric. In the following we 
can set Aq = 1 as this does not affect the results. There are several options here , 
because wc need to pass from the Einstein to the string metric and this depends on 
whether is the string dilaton or not. We parametrize this dependence by 

= bE eV^"^ (C.ll) 

with k = 1 corresponding to the string dilaton and A; = to a completely different 
scalar. 

We obtain _ 

bs^(-^ , d = a-^J^^L= (C.12) 

The turning point Ws of the dragging solution is given by 

f{w,) = v' ^ = Wo{l - v'f^ ^ ^{1 - v'f^ (C.13) 

Prom now we use the string frame scale factor bg that we still call b. 
We can compute 

^ = -A = 4.r2(i-«ri , ^^—d-v')-^ (C.14) 

and 



(C.15) 
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Note that as ac > 0, the last square root never vanishes. 

The appropriate couphng to the dilaton is the one that makes 4a = c, and then 



k 



v/2(p - 1) 2{p - If 5^ - 7{p - 1)S^ - 2{p - 4) 
4 2 - (p - 1)52 

We now move to study the equations for the two point functions 



(C.16) 



a. 



du, 
with 



rWif-v'Wf-c') (5x11) 



Z2 



+ 



6X^ = (C.17) 
5X11 = (C.18) 



(C.19) 



Define 



^ d^^—d, , f{w)-v'^{l-v'){l-x^) (C.20) 

Ws Wo Ws 



h\w)f - = 



4a 



1 - (1 - i;2)x^ - v'^x^^ 



X 



4a 



X'- 



1 - (1 - v'^)x'' - v^x^" 



(C.21) 



(C.22) 



We then have 



x-^V(l - - (1 - - v^x^'')d^ (5X^) 

^2 



+ 



with 



^(1 - X^){\ - (1 - - y22,4a) 



Ws w c(l — 2c a; 



=5X^ = (C.23) 



47r T 



(C.24) 



For 5 = 0, p = 4, /c = the dependence on v cancels and it reduces to the AdS^ case. 
When 4a = c, = 1 — and the equation simplifies to 



X 



Q2 
(1 - X'^) 



(5X^) = (C.25) 



-46- 



References 



[1] J. Adams et al. [STAR Collaboration], "Experimental and theoretical challenges in 
the search for the quark gluon plasma: The STAR collaboration's critical assessment 
of the evidence from RHIC collisions, " Nucl. Phys. A 757 (2005) 102 



ArXiv:nucl-ex/0501009 



B. B. Back et al., "The PHOBOS perspective on discoveries at RHIC," Nucl. Phys. 
A 757 (2005) 28 ||ArXiv:nucl-ex/0410022| |; 

I. Arsene et al. [BRAHMS Collaboration], "Quark gluon plasma and color glass 
condensate at RHIC? The perspective from the BRAHMS experiment, " Nucl. Phys. 
A 757 (2005) 1 ||ArXiv:nucl-ex/0410020|]; 



K. Adcox et al. [PHENIX Collaboration], "Formation of dense partonic matter in 
relativistic nucleus nucleus collisions at RHIC: Experimental evaluation by the 
PHENIX collaboration, " Nucl. Phys. A 757 (2005) 184 ||ArXiv:nucl-ex/0410003|| . 



[2] D. T. Son and A. O. Starinets, "Viscosity, Black Holes, and Quantum Field 



Theory, " Ann. Rev. Nucl. Part. Sci. 57 (2007) 95 | ArXiv:0704.0240| |hep-th 



[3] J. Casalderrey-Solana and C. A. Salgado, "Introductory lectures on jet quenching in 
heavy ion collisions, " Acta Phys. Polon. B 38 (2007) 3731 1| ArXiv:0712.3443| |hep-ph1 . 



[4] E. lancu, "Partons and jets in a strongly- coupled plasma from AdS/CFT," Acta 
Phys. Polon. B 39 (2008) 3213 || ArXiv:0812.0500| Ihep-phj . 



[5] S. S. Gubser, S. S. Pufu, F. D. Rocha and A. Yarom, "Energy loss in a strongly 
coupled thermal medium and the gauge-string duality, " |[ArXiv:0902.404r hep-th]| . 



[6] U. Gursoy, E. Kiritsis, L. Mazzanti, G. Michalogiorgakis and F. Nitti, "Improved 
Holographic QCD, " || ArXiv: 1006 . 546 1 1 1 hep-thj . 



[7] C. P. Herzog, A. Karch, P. Kovtun, C. Kozcaz and L. G. Yaffe, "Energy loss of a 
heavy quark moving through N = 4 supersymmetric Yang-Mills plasma, " JHEP 0607 



(2006) 013; ||ArXiv:hep-th/0605158 



H. Liu, K. Rajagopal and U. A. Wiedemann, "Calculating the jet-quenching 
parameter from AdS/CFT, " Phys. Rev. Lett. 97 (2006) 182301 
ArXiv:hep-ph/0605178|| ; 



"Wilson loops in heavy ion collisions and their calculation in AdS/CFT, " JHEP 
0703 (2007) 066 ||ArXiv:hep-ph/0612168|| . 

[9] S. S. Gubser, "Drag force in AdS/CFT," Phys. Rev. D 74 (2006) 126005 
II ArXiv:hep-th/0605 182]| 



[10] H. Liu, K. Rajagopal and Y. Shi, "Robustness and Infrared Sensitivity of Various 
Observables in the Application of AdS/CFT to Heavy Ion Collisions," JHEP 0808 



(2008) 048 ArXiv:0803.3214 hep-ph 



-47- 



[11] 



[121 



U. Gursoy, E. Kiritsis, G. Michalogiorgakis and F. Nitti, "Thermal Transport and 
Drag Force in Improved Holographic QCD, " JHEP 0912 (2009) 056 



ArXiv:0906.1890][hep-ph" 



J. Casalderrey-Solana and D. Teaney, "Heavy quark diffusion in strongly coupled N 
= 4 Yang Mills, " Phys. Rev. D 74 (2006) 085012 ||ArXiv:hep-ph/0605199|| . 



C. P. Herzog and D. T. Son, "Schwinger-Keldysh propagators from AdS/CFT 
correspondence, " JHEP 0303, 046 (2003) || ArXiv:hep-th/0212072 1 . 



[13] 

[14] S. S. Gubser, "Momentum fluctuations of heavy quarks in the gauge-string duality," 
Nucl. Phys. B 790 (2008) 175 |ArXiv:hep-th/0612143] . 

[15] J. Casalderrey-Solana and D. Teaney, "Transverse momentum broadening of a fast 
quark in a N = 4 Yang Mills plasma, " JHEP 0704 (2007) 039 



ArXiv:hep-tli/0701123] 



[16] J. de Boer, V. E. Hubeny, M. Rangamani and M. Sliigemori, "Brownian motion in 
AdS/CFT," JHEP 0907 (2009) 094 |[ArXiv:0812.5112] [li^^^ . 

[17] D. T. Son and D. Teaney, "Thermal Noise and Stochastic Strings in AdS/CFT," 
JHEP 0907, 021 (2009) || ArXiv:0901.2338| [liep-th| . 

[18] G. C. Giecold, E. lancu and A. H. Mueller, "Stochastic trailing string and Langevin 
[19] 



dynamics from AdS/CFT, " JHEP 0907, 033 (2009) || ArXiv:0903.1840| [hep-tti 



E. Caceres, M. Chernicoff, A. Guijosa and J. F. Pedraza, "Quantum Fluctuations 
and the Unruh Effect in Strongly- Coupled Conformal Field Theories," JHEP 1006, 



078 (2010) II ArXiv:1003. 53321 |hep-th 



[20] 
[21] 

[22] 
[23] 

[24] 



C. Hoyos-Badajoz, "Drag and jet quenching of heavy quarks in a strongly coupled 
N =2* plasma," JHEP 0909, 068 (2009) |[ArXiv:arXiv:0907.5036][hep-tliJ . 



U. Gursoy, E. Kiritsis, L. Mazzanti and F. Nitti, "Langevin diffusion of heavy quarks 
in non-conformal holographic backgrounds, " JHEP 1012, 088 (2010) 



[ArXiv: 1006.3261] [hep-th] 



Y. Akiba [PHENIX Collaboration], "Probing the properties of dense partonic matter 
at RHIC," Nucl. Phys. A 774 (2006) 403 [arXiv:nucl-ex/0510008]. 

S. S. Adler et al. [PHENIX Collaboration], "Nuclear modification of electron spectra 
and implications for heavy quark energy loss in Au -h Au collisions at s(NN)**(l/2) 
= 200-GeV," Phys. Rev. Lett. 96 (2006) 032301 || ArXiv: nucl-ex/05 10047 1| . 

B. I. Abelev et al. [STAR Collaboration], "Transverse momentum and centrality 
dependence of high-px non-photonic electron suppression in Au-hAu collisions at 
y^sj^ = 200 GeV," Phys. Rev. Lett. 98 (2007) 192301 ||ArXiv:nucl-ex/0607012 



-48 - 



A. Adare et al. [PHENIX Collaboration], "Energy Loss and Flow of Heavy Quarks in 
Au+Au Collisions at ^/I]^ = 200 GeV," Phys. Rev. Lett. 98 (2007) 172301 



ArXiv:nucl-ex/0611018 



T. C. Awes [PHENIX Collaboration], "Highlights from PHENIX - U," J. Phys. G 35 
(2008) 104007 II ArXiv: 0805 ■ 1636 1 1 nucl-ex] . 



N. Armesto, M. Cacciari, A. Dainese, C. A. Salgado and U. A. Wiedemann, "How 
sensitive are high-p(T) electron spectra at RHIC to heavy quark energy loss?," Phys. 
Lett. B 637 (2006) 362 || ArXiv:hep-ph /05 1 1 257 1 . 



R. Rapp and H. van Hees, "Heavy Quarks in the Quark- Gluon Plasma," 



ArXiv:0903.1096 hep-ph 



Y. Akamatsu, T. Hatsuda and T. Hirano, "Heavy Quark Diffusion with Relativistic 



Langevin Dynamics in the Quark- Gluon Fluid," [ArXiv: 0809. 1499] [hep-ph 



F. Debbasch, K. Malhck and J.P. Rivet, J. of Stat. Phys. 88, 945 (1997). 
F. Debbasch and J.P. Rivet, J. of Stat. Phys. 90, 1179 (1998). 
C. Chevaher and F. Debbasch, J. Math. Phys. 49, 043303 (2008). 

J. Dunkel, P. Hanggi, Relativistic Brownian Motion, Phys. Rep. 471 (2009) 1. 



ArXiv: 08 1 2 . 1996 cond-mat .stat-mech 



S. S. Gubser, "Comparing the drag force on heavy quarks in N = 4 super-Yang-Mills 



theory and QCD, " Phys. Rev. D 76 (2007) 126003 ||ArXiv:hep-th/0611272 



U. Gursoy and E. Kiritsis, "Exploring improved holographic theories for QCD: Part 



I, " JHEP 0802 (2008) 032 ArXiv:0707.1324 hep-th 



U. Gursoy, E. Kiritsis and F. Nitti, "Exploring improved holographic theories for 



QCD: Part H, " JHEP 0802 (2008) 019 ||ArXiv:0707.1349| [hep-th 



S. S. Gubser and A. Nellore, "Mimicking the QCD equation of state with a dual black 



hole, " Phys. Rev. D 78 (2008) 086007 ArXiv:0804.0434 hep-th 



U. Gursoy, E. Kiritsis, L. Mazzanti and F. Nitti, "Deconfinement and Gluon Plasma 
Dynamics in Improved Holographic QCD," Phys. Rev. Lett. 101, 181601 (2008) 



ArXiv:0804.0899 hep-th 



U. Gursoy, E. Kiritsis, L. Mazzanti and F. Nitti, "Holography and Thermodynamics 



of 5D Dilaton- gravity, " JHEP 0905 (2009) 033 |[ArXiv:0812.0792 |[hep-tli 



E. Kiritsis, "Dissecting the string theory dual of QCD, " Fortsch. Phys. 57 (2009) 396 



ArXiv:0901.T772] hep-th 



U. Gursoy, E. Kiritsis, L. Mazzanti and F. Nitti, "Improved Holographic QCD at 
finite Temperature: comparison with data, " Nucl. Phys. B 820 (2009) 148 



ArXiv:0903. 28591 [hep-th 



- 49 - 



M. Panero, "Thermodynamics of the QCD plasma and the large-N limit, " Phys. Rev. 
Lett. 103 (2009) 232001 || ArXiv:0907.3719| |hep-lat] 



[40] 
[41] 

[42] A. Karch and A. O'Bannon, "Metallic AdS/CFT," JHEP 0709 (2007) 024 



E. Kiritsis, "Supergravity, D-brane probes and thermal super Yang-Mills: A 
comparison," JHEP 9910 (1999) 010 ||ArXiv:hep-th/9906"206i 



ArXiv:0705.3870] [hep-th 



[43] 
[44] 
[45] 



S. Janiszewski, A. Karch, "Moving Defects in AdS/CFT, " ||ArXiv:1106.4010||hep-"th1 . 
J. Casalderrey-Solana, D. Fernandez and D. Mateos, "A New Mechanism of Quark 



Energy Loss, " Phys. Rev. Lett. 104 (2010) 172301 |ArXiv:0912.3717| |hep-ph 



L. F. Cughandolo, J. Kurchan and L. Pehti, "Energy flow, partial equilibration, and 
effective temperatures in systems with slow dynamics," Phys. Rev. E 55 (1997) 3898. 

D. Cubero, J. Casado-Pascual, J. Dunkel, P. Talkner, P. Hanggi, "Thermal 
equilibrium and statistical thermometers in special relativity" Phys. Rev. Lett. 



99:170601 (2007); [[ArXiv:0705.3328] [cond-mat.stat-mech 



[46] N. Iqbal and H. Liu, "Universality of the hydrodynamic limit in AdS/CFT and the 
[47] 



membrane paradigm, " Phys. Rev. D 79 (2009) 025023 [ArXiv:0809.3808J [hep-th 



D. Teaney, "Finite temperature spectral densities of momentum and R- charge 
correlators in N = 4 Yang Mills theory, " Phys. Rev. D 74 (2006) 045025 



ArXiv:hep-ph/0602044 



[48] 

[49] 

[50] 
[51] 
[52] 

[53] 
[54] 



J. Casalderrey-Solana and C. A. Salgado, "Introductory lectures on jet-quenching in 
heavy ion collisions, " Acta Phys. Polon. B 38, 3731 (2007) 



ArXiv:0801.2173][hep-ph 



R. P. Feynman and F. L. . Vernon, "The theory of a general quantum system 
interacting with a linear dissipative system," Annals Phys. 24, 118 (1963) [Annals 
Phys. 281, 547 (2000)]. 

H. Kleinert, 'T*ATH INTEGRALS in Quantum Mechanics, Statistics, Polymer 
Physics, and Financial Markets , " ^PIRES entry World Scientific, Singapore, 2004. 



S. S. Gubser, "Comparing the drag force on heavy quarks in N = 4 super- Yang- Mills 
theory and QCD, " Phys. Rev. D 76 (2007) 126003 ||ArXiv:hep-th/0611272|| . 



D. T. Son and A. O. Starinets, "Minkowski- space correlators in AdS/CFT 
correspondence: Recipe and applications, " JHEP 0209, 042 (2002) 



ArXiv:hep-th/0205051 



S. S. Gubser, S. S. Pufu and F. D. Rocha, "Bulk viscosity of strongly coupled plasmas 
with holographic duals, " JHEP 0808, 085 (2008) || ArXiv:0806.0407[ Ihep-th] . 



K. B. Fadafan, H. Liu, K. Rajagopal and U. A. Wiedemann, "Stirring Strongly 
Coupled Plasma, " Eur. Phys. J. C 61 (2009) 553 ||ArXiv:0809.2869/|hep-p"h|' 



- 50 - 



[55] P. M. Hohler, M. A. Stephanov, "Bulk spectral function sum rule in QCD-like 
theories with a holographic dual," JHEP 1106, 130 (2011). 



ArXiv:1103.0977][hep-ph 



[56] S. Caron-Huot, P. M. Chesler and D. Teaney, "Fluctuation, dissipation, and 
thermalization in non- equilibrium AdS^ black hole geometries, " 



ArXiv:1102.1073||hep-th 



[57] C. Charmousis, B. Gouteraux, B. S. Kim, E. Kiritsis and R. Meyer, "Effective 

Holographic Theories for low-temperature condensed matter systems," JHEP 1011 



(2010) 151 |[ArXiv: 1005.4690] [hep-th] 



[58] I. Kanitscheider, K. Skenderis and M. Taylor, "Precision holography for 
non-conformal branes, " JHEP 0809 (2008) 094 || ArXiv:0807.3324| |hep-th| ; 
I. Kanitscheider and K. Skenderis, "Universal hydrodynamics of non-conformal 



branes, " JHEP 0904 (2009) 062 || ArXiv:0901.1487| Ihep-tlT 



[59] B. Gouteraux and E. Kiritsis, "Generalized Holographic Quantum Criticality at 
Finite Density, " |[ArXiv:1107.2116] [hep-th]] . 

[60] F. Bigazzi, A. L. Cotrone, J. Mas, A. Paredes, A. V. Ramallo and J. Tarrio, "D3-D7 
Quark-Gluon Plasmas," JHEP 0911, 117 (2009) |[ArXiv: 0909. 2865] [hep-th] . 

[61] C. Nunez, A. Paredes and A. V. Ramallo, "Unquenched flavor in the gauge/gravity 
correspondence, " |[ArXiv: 1002. 1088 hep-th]| . 

[62] G. Bertoldi, F. Bigazzi, A. L. Cotrone and J. D. Edelstein, "Holography and 
Unquenched Quark-Gluon Plasmas," Phys. Rev. D 76, 065007 (2007) 



ArXiv:hep-th/0702225 



[63] F. Karsch, D. Kharzeev and K. Tuchin, "Universal properties of bulk viscosity near 
the QGD phase transition," Phys. Lett. B 663 (2008) 217 ||ArXiv:0711.0914| |hep-ph|| . 

[64] Y. Kinar, E. Schreiber, J. Sonnenschein, "Q anti-Q potential from strings in curved 
space-time: Classical results, " Nucl. Phys. B566 (2000) 103-125. 



ArXiv:hep-th/9811192 



[65] E. Kiritsis, L. Mazzanti and F. Nitti, to appear. 



- 51 - 



